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Preface 


The  purpose  of  this  study  was  to  examine  the  mode 
eigenvalues,  far  field  beam  quality,  and  far  field  beam 
steering  of  decentered,  unstable  resonators.  Because  of 
the  complex  nature  of  mode  modeling,  much  of  the  anlaysis 
performed  was  numerical  instead  of  theoretical.  It  is 
hoped  that  the  lack  of  explicit  formulae  and  mathematical 
developments  will  not  detract  from  the  value  of  the 
numerical  analyses. 

As  with  any  major  undertaking,  this  study  would  not 
have  been  possible  without  the  support,  aid,  and  encourage¬ 
ment  of  many  people.  Much  gratitude  and  appreciation  is 
due  Lt  Col  John  Erkkila,  my  advisor,  for  his  constant  guid¬ 
ance  and  understanding,  especially  during  the  final  phases 
of  this  project.  My  thanks  are  due  Capt  Mark  Rogers  for 
the  many  useful  discussions  of  resonators  and  his  help  with 
the  computer  (a  monster  by  any  measure!).  Sharon  Gabriel, 
my  typist,  deserves  much  appreciation  for  her  excellent 
work.  And  finally,  my  sincerest  appreciation  goes  to  my 
family  and  friends  for  their  support  and  words  of  encourage¬ 
ment  throughout  this  endeavor. 


Steven  M.  Rinaldi 
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Abstract 

r  » 

The  mode  eigenvalues,  far  field  integrated  intensity, 
and  far  field  beam  steering  angles  of  unstable,  decentered 
strip  resonators  were  studied .  The  resonators  examined  had 
magnifications  of  2.0  and  equivalent  Fresnel  numbers  in  the 

-  .*  l  • 

range  9.3  <  N  <9.9.  The  resonator  modes  were  calculated 
by  the  asymptotic  method  of  Horwitz. 

Two  equivalent  Fresnel  numbers  for  the  decentered  reson¬ 
ators  were  defined.  The  fundamental  and  second-order  mode 
eigenvalues  exhibited  periodicities  in  the  equivalent  Fresnel 
numbers.  The  mode  separation  was  observed  to  be  a  function 
of  the  amount  of  decenter  and  the  two  equivalent  Fresnel  num¬ 
bers.  The  cusps  of  the  first  two  eigenvalues  were  cyclic  in 

N  . 
eq 

The  far  field  integrated  intensity  was  computed  for  spot 
sizes  of  one,  two,  and  three  Airy  disks.  The  percentage  of 
total  power  deposited  in  a  given  spot  size  increased  as  the 
decenter  increased.  Beam  quality  instabilities  were  observed 
in  all  modes. 

The  beam  steering  angles  of  the  first  four  modes  were 
calculated.  The  angles  fluctuated  about  the  optic  axis  as  the 
decenter  was  increased.  The  fundamental  mode  had  significantly 
lower  beam  steering  than  the  higher-order  modes. 
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PERFORMANCE  ANALYSIS  OF 


DECENTERED  UNSTABLE  RESONATORS 

I.  Introduction 


Background 

Optical  cavities  may  be  categorized  in  two  general 
classes:  the  stable  and  the  unstable  resonators  Stable 
resonators  are  characterized  by  well-defined  mod'  The 
mode  volumes  of  such  cavities  are  generally  quit  .•'■all. 
Output  coupling  of  the  modes  is  usually  accomplifc...~d  by 
transmission  through  a  mirror  or  other  optical  element. 
Consequently,  the  output  power  levels  achievable  in  stable 
resonator  designs  are  limited  to  relatively  low  values. 
Unstable  resonators  are  characterized  by  large  mode  volumes. 
The  modes  are  outcoupled  via  diffraction  around  one  or  both 
of  the  resonator  mirrors  The  output  power  levels  of 
unstable  resonator  lasers  are  not  limited  to  low  values; 
hence,  unstable  resonators  are  used  in  high  power  applica¬ 
tions  . 

The  modes  in  a  stable  resonator  are  well-defined  and 
easily  calculated.  The  modes  are  generally  described  in 
terms  of  the  Hermite-Gaussian  or  Laguerre-Gaussian 
functions  (Ref  1:1324).  The  modes  of  unstable  resonators 
are  not  so  simply  described.  They  must  be  computed  by 


using  one  of  several  different  numerical  techniques.  One 
solution  technique  of  particular  interest  is  the  asymptotic 
method  (Ref  2) ,  which  was  used  to  determine  the  modes  in 
this  study. 

The  geometry  of  a  general  resonator  is  depicted  in 
Figure  1-1.  M  is  the  geometric  magnification,  a^  and 
a 2  are  the  linear  half-widths  of  mirrors  and  M2  , 

respectively,  R^  and  R^  are  the  radii  of  curvature  of 
the  mirrors,  and  L  is  the  axial  length  of  the  resonator. 
The  well-known  g  parameters  are  given  by 


R. 


If  0  <  g-j^  s  1  #  fche  resonator  is  stable.  Otherwise, 

the  resonator  is  classified  as  unstatle. 

The  equivalent  Fresnel  number  is  defined  by 


2X 'Lgx 


(1 


X'  is  the  radiation  wavelength.  This  parameter  is  related 
to  the  number  of  Fresnel  zones  intercepted  by  mirror 


The  geometry  of  misaligned  and  decentered  resonators 
is  depicted  in  Figure  1-2.  It  can  easily  be  shown  (Ref  3: 
19-23)  that  a  misalignment  of  mirror  M2  is  equivalent  to 
a  decenter  of  that  mirror: 
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Figure  1-1.  Geometry  of  an  Unstable  Laser  Resonator 


a2  0  M  +  1 


where  6  is  the  fractional  offset  of  and  0  is  the 

angle  through  which  is  tilted.  As  a  result,  off-axis 

resonators  can  be  analyzed  as  misaligned  resonators,  and 
vice-versa . 

The  confocal  resonator  is  a  special  type  of  cavity. 
The  g  parameters  of  this  resonator  must  satisfy  the 
relationship : 


Ml. ALIGNED  RESONATOR 


DECENTERED  RESONATOR 


Figure  1-2.  Geometries  of  Misaligned 
and  Decentered  Resonators 


gl  +  g2  =  2glg2 


(1.4) 


The  output  beam  of  a  confocal  resonator  is  collimated.  In 
the  geometric  optics  limit,  the  beam  is  a  plane  wave.  This 
is  a  limiting  value  of  the  spherical  wave  mode  (in  the 
geometric  optics  limit)  of  a  general  unstable  resonator. 

The  Asymptotic  Method 

This  section  discusses  the  asymptotic  method  of 
solving  for  bare,  strip  resonator  modes.  As  this  solution 
technique  is  derived  in  detail  elsewhere  (Refs  2,  3,  4,  5), 
the  following  development  is  simply  a  brief  outline. 

The  generalized  integral  equation  describing  the 
resonator  modes  is 


b 

Au  (x)  =  /  K  (x  ,y)  u  (y)  dy  (1.5) 

a 


where  u(x)  is  the  mode,  A  is  the  mode  eigenvalue,  and 
K(x,y)  is  the  kernel  of  the  integral.  For  a  decentered 
resonator,  Eq  (1.5)  becomes 


Ag  (x) 


exp [ -it (y-x/M) 2 ] dy 


(1.6) 
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where 


inN  x2 

g  (x)  =  e  ec^  u(x) 


t  =  tiMF 


irMa2 

2A 


(1.7) 


(1.8) 


In  Eq  (1.6),  the  half-width  a2  of  mirror  M2  has  been 
normalized  to  unity.  The  integration  is  thus  taken  over 
the  surface  of  M2  ,  with  the  decenter  accounted  for. 
y  is  a  dummy  variable  of  integration. 

g (x)  can  be  approximated  as  a  unit  amplitude  spherical 
wave  plus  a  finite  series  of  higher-order  edge-diffracted 
waves.  This  expansion  is  given  by 


N 


g(x)  =  1+1  (anFn(x)  +  bnGn(x)] 

n=l 


(1.9) 


where 


Fn(x) 


n-l  exp[-it(l-x/M)  /Mn_] ] 


4iiTt 


(1  -  x/M") 


(1.10) 


m 


Gn(x)  = 


n.  2 


n_!  exp  [-it  (1+x/M  ) 


4iirt 


(1  +  x/Mn) 


(1.11) 


M 


n 


n 

l  M 

k=0 


-2k 
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Equations  (1.9),  (1.10),  and  P.ll)  are  substituted 
into  Eq  (1.6).  The  resultant  integral  equation  is  evaluated 


by  the  method  of  stationary  phase.  The  final  result  is 


N 

X { 1  +  l  [a  F  (x)  +  b  G  (x)  ]  }  =  1  +  F.  (x)  +  G.  (x) 

_  n  n  nn  1  1 


N 

+  F  [a  F  ,  (x) 
1  n  n+1 

n=l 


bnGn+l(x)l 


+  F1(x) 


N 


l 

n=l 


[a  F  ( B)  +  b  G  ( B)  ] 
n  n  n  n 


N 

+  G.  (x)  l 
1  n=l 


[a  F  (a)  +  b  G  (a)  ] 
n  n  n  n 


(1.13) 


where 

a  =  -1+6  (1.14) 

6  =  1+6  (1.15) 

As  shown  in  Chapter  II,  Eq  (1.13)  can  be  reduced  to 

a  polynomial  expression  in  X  with  known  coefficients.  The 

equation  can  be  solved  numerically  for  the  eigenvalues.  The 

coefficients  a  and  b  are  then  easily  computed, 
n  n 

Finally,  the  mode  u(x)  is  calculated  from  Eqs  (1.9)  and 
(1.7)  . 
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Objectives 

The  objectives  of  this  study  are  to  examine  the 
behavior  of  the  eigenvalues  and  the  far  field  modes  as 
functions  of  the  decenter  parameter  6  .  Specifically, 

the  objectives  are: 

(1)  To  determine  how  the  eigenvalues  evolve  as  6 

is  increased  from  zero.  The  relationships 

between  the  cusping  nature  of  the  first  two 

eigenvalues  and  the  equivalent  Fresnel  numbers 

N  r  and  N  ..  will  be  explored  (see 
eq,L  eq,U  c 

Chapter  II) . 

(2)  To  determine  how  the  far  field  integrated 
intensity  changes  as  a  function  of  6  for 
the  first  four  modes. 

(3)  To  explore  the  relationship  between  the  far 
field  beam  steering  angles  and  6  for  the 
first  four  modes. 

(4)  To  develop  a  simple  set  of  design  criteria 
for  unstable  resonators  based  on  the  observa¬ 
tions  listed  above. 

Assumptions  and  Limitations 

The  assumptions  and  limitations  of  Reference  3  apply 


to  this  study,  as  the  computer  model  developed  in  that 
work  was  used  to  predict  the  resonator  modes.  The  following 
additional  assumptions  and  limitations  are  made: 


(1)  The  Fraunhofer  approximation  is  sufficient  to 
calculate  the  far  field  mode  patterns.  This 
approximation  is  valid  if  the  far  field  patterns 
are  calculated  at  distances  far  from  the  resonator 
output  aperture.  Since  this  separation  was  taken 
to  be  infinite,  this  approximation  is  valid. 

(2)  The  computer  code  used  to  predict  the  resonator 
modes  is  valid  for  decenters  of  0.0  <  o  <  0.9. 
Weiner  asserts  that  this  is  essentially  true 
(Ref  6:1831),  based  on  comparisons  of  modes 
calculated  by  the  asymptotic  method  and  the 
power  method  (Ref  7) . 

(3)  Only  bare  strip  resonators  with  magnification 

M  =  2  and  equivalent  Fresnel  numbers  in  the 

range  9.3  <  N  <  9.9  will  be  studied.  The 
eq 

results  of  Appendix  A,  however,  are  valid  for 

general  M  and  N  values. 

eq 

Throughout  this  study,  reference  will  be  made  to  the 
"geometric  mode."  This  is  the  geometric  mode  of  a  confocal 
resonator  -  a  uniform  plane  wave.  Upon  emergence  from  the 
resonator,  the  mode  will  have  a  uniform  "annular"  profile. 

Organization 

In  Chapter  II,  a  general  theory  of  mode  eigenvalues 


is  presented.  The  eigenvalue  polynomial  equation  is  devel¬ 
oped,  and  theoretical  predictions  about  eigenvalues  of 


decentered  resonators  are  made.  The  far  field  integrated 
intensity  and  beam  steering  are  discussed  in  Chapter  III, 
based  on  a  Fourier  optics  treatment  of  the  geometric  mode. 

In  Chapter  IV,  the  major  sources  of  numerical  error  that 
exist  in  the  computer  code  that  determines  the  far  field 
modes,  integrated  intensity,  and  beam  steering  are  evaluated. 
The  study  results  are  presented  in  Chapter  V.  The  general 
unstable  resonator  design  criteria  are  presented  in  Chapter 
VI,  and  the  overall  conclusions  and  recommendations  are 
relegated  to  Chapter  VII. 
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II .  Theory  of  Mode  Eigenvalues 

A  general  theory  of  the  mode  eigenvalues  is  presented. 

The  polynomial  equation  used  to  compute  the  eigenvalues 
in  the  asymptotic  approximation  is  derived  for  completeness, 
even  though  its  development  exists  elsewhere  (Refs  3:28-30; 
•4-:24-29)  .  The  concept  of  two  equivalent  Fresnel  numbers 
for  decentered  resonators  is  developed,  and  its  relationship 
to  the  mode  eigenvalues  is  discussed. 

Background  Theory  of  Eigenvalues  -  Aligned  Resonators 

The  general  form  of  the  integral  equation  describing 
the  round  trip  propagation  through  a  resonator  is 

b 

Xu(x)  =  /  K (x,y) u (y) dy  (2.1) 

a 

This  is  an  eigenvalue  problem.  X  is  the  eigenvalue  and 
the  integration  is  the  operator. 

After  a  round  trip  through  the  resonator  under  steady 
state  conditions,  the  mode  must  be  essentially  unchanged  in 
form.  After  the  propagation,  the  mode  u/(x)  must  be 
given  by 


u' (x)  =  Xu(x) 


(2.2) 


At  every  point  on  the  wavefront,  the  amplitude  is  scaled 


11 


by  the  magnitude  of  the  eigenvalue  and  the  phase  is 
shifted  by  the  phase  of  the  eigenvalue. 


The  fraction  of  energy  coupled  out  of  the  resonator 
is  related  to  the  eigenvalue.  For  strip  resonators. 


1  i  ‘ 

(Outcoupled  Energy  Fraction) ^  =  1  -  — — —  (2.3) 


a.  t-  • 

where  i  refers  to  the  i  mode  and  M  is  the  geometric 
magnification  of  the  resonator. 

The  symmetric  mode  eigenvalues  exhibit  a  periodicity 
when  plotted  against  the  equivalent  Fresnel  number  of  the 
resonator.  Notably,  the  separation  between  the  magnitudes 
of  the  first  and  second  symmetric  mode  eigenvalues  has 
maxima  at  approximately 


Neq  n  +  8  '  n-0,1,2, . . .  (2.4) 

and  minima  at  approximately 

Neq  =  n  +  I  ,  n=0, 1, 2 , . . .  (2.5) 


This  behavior  is  apparent  in  Figure  2-1.  This  plot  displays 
the  magnitudes  of  the  first  seven  (symmetric  and  anti¬ 
symmetric)  mode  eigenvalues.  A  periodic  nature  of  the 
eigenvalues  is  rather  clear. 
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Below  some  critical  equivalent  Fresnel  number, 

Neq(crit)  '  t*le  ^rst  an<^  second  symmetric  eigenvalues 

display  either  crossing  or  cusping  at  the  mode  separation 

minima  points.  Above  Neq(crj_t)  '  crossin9  of  the  first 

two  eigenvalues  ceases.  N  ,  . ..  is  given  by  (Ref  9:4149): 

eq(cnt) 


Neq (crit) 


11.5 
(in  M)  3 


(2.6) 


Above  Neq(crit)  '  the  overall  periodicity  of  the  eigen¬ 
values  still  exists.  The  amplitude  of  the  fluctuations 
of  the  magnitude  of  decrease* as  N  increases  (Ref  2: 

1536)  . 

The  periodic  fluctuations  of  |X^|  may  be  better 
understood  if  the  following  argument  is  advanced.  The 
equivalent  Fresnel  number  is  equal  to  the  number  of  half 
wavelengths  between  the  edge  of  the  output  mirror  and  the 
closest  point  on  the  geometric  wave  inside  the  resonator 
when  that  wave  just  touches  the  center  of  the  mirror 
(Ref  10:360).  This  is  depicted  in  Figure  2-2.  Of  all  the 
waves  that  are  scattered  from  the  edge  of  the  mirror,  that 
which  is  propagated  back  into  the  resonator  along  the  ray 
direction  of  the  outgoing  geometric  wave  is  the  most 
important.  It  is  focused  back  along  the  optic  axis,  where 
it  can  interfere  with  the  resonator  mode  (Ref  11:263).  As 
N  is  changed  by  unity,  the  phase  of  this  re-entrant  ray 
is  changed  by  2ir  radians.  Consequently,  there  should  be 
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Figure  2-2.  Interpretation  of  N  (Ref  10:360) 

eg 

cyclical  constructive  and  destructive  interference  along 
the  axis,  with  a  period  of  unity  in  . 

Considering  just  the  fundamental  mode,  when  the  re¬ 
entrant  ray  interferes  destructively  with  the  mode,  the 
field  on  the  resonator  axis  should  be  relatively  diminished 
compared  to  the  field  near  the  edges  of  the  mode.  If 
the  re-entrant  ray  interferes  constructively  with  the  mode, 
the  field  on  the  optic  axis  should  be  relatively  more 
intense  than  the  field  near  the  edges  of  the  mode.  Hence, 

a  resonator  with  N  such  that  destructive  interference 

eq 

occurs  should  outcouple  relatively  more  power  than  a 

resonator  with  N  such  that  constructive  interference 

eq 

occurs . 

Since  the  outcoupled  energy  of  a  mode  is  related  to 
| X |  as  in  Eq  (2.3),  | A |  should  also  be  related  to  the 
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constructive  and  destructive  interference.  When  the  re¬ 
entrant  ray  interferes  constructively,  |X|  should  be 
near  its  maximum.  When  destructive  interference  occurs, 

| X |  should  be  near  its  minimum.  Hence,  for  |A|  near 
its  maximum,  the  mode  should  be  relatively  built  up  on 
the  resonator  axis  and  depressed  near  its  edges.  When 
| X |  is  near  a  minimum  value,  the  mode  should  be  depressed 
on  the  optic  axis  and  built  up  near  its  edges. 

An  examination  of  the  fundamental  mode  for  5  <  N  <21 

eq 

confirms  the  above.  1^-jJ  reaches  its  peak  value  when 
3 

N  =  n  +  g-  (n=0 , 1 , 2  ,  . .  . )  ,  and  its  minimum  value  when 

N  =  n  +  Z-  .  The  fundamental  mode  does  have  a  depressed 

eq  8 

,  ~  7 

intensity  on  axis  when  N  =  n  +  ^  and  relative  intensity 
2  eq  8  2 

3 

peaks  on  axis  when  =  n  +  g  .  Figures  2-3  and  2-4  show 

representative  plots  of  this  behavior. 


The  stationary  phase  approximation  to  the  integral 
equation  results  in  the  expression 


N 

A { 1  +  l  [anFn(x)  +  b  G  (x)  ]  }  =  1  +  F^x)  +  G1(x) 
n=l 
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With  the  manipulations  detailed  below,  Eq  (2.7)  can 
be  reduced  to  a  polynomial  in  X  of  order  2N+1  .  The 

polynomial  can  be  numerically  solved  for  the  2N+1  eigen¬ 
values. 

The  initial  step  is  to  equate  the  constants  and  the 

coefficients  of  F  and  G  on  both  sides  of  Eq  (2.7). 

n  n 

Equating  the  constants  yields 


X 


1  + 


aNFN+l  + 


bNGN+l 


(2.8) 


Fn+i  and  are  essentially  constant,  as  may  be  seen 

in  Eqs  (1.10)  and  (1.11) . 

Equating  the  coefficients  of  F^x)  and  (x) 

gives 

N 

Xai  =  1  +  l  [anFn(6)  +  bnGn(B)]  (2.9a) 

n=  1 

N 

Xbx  =  1  +  I  [anFn(cO  +  bnGn(a)]  (2.9b) 

n=l  '  “ 

Finally,  equating  the  coefficients  of  F  (x)  and  Gn(x)  / 
n^l  ,  yields 


Xa 


n+1 


a 

n 


(2.10a) 


Xb 


n+1 


b 


n 


(2.10b) 
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From  Eq  (2.10a), 


_n  _  _1_ 

n+l  A  .  n 


(2.11) 


Changing  the  subscripts  slightly,  and  using  Eq  (2.11)  twice 
yields 


a  A 
n 


n-1 


,  N-l 


N-1 


(2.12) 


=  a  A 
n 


n-N 


(2.13) 


Similarly, 


b  =  b  A 
N  n 


n-N 


(2.14) 


Using  Eqs  (2.13)  and  (2.14),  Eq  (2.9)  reduces  to 


V 


N 


N 


1  +  l,  jN'"  IVn1'1  +  bnV8,l 

n=l 


(2.15a) 


bNX 


N 


N 


1  *  l  XN'n  [aNFn(0  +  bsGn(c)] 
n=l 


(2.15b) 


Following  Reference  3,  define 
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'  =  )  A"  "F  ot 

a  L ,  n 

n=l 


I  AN'nF  (B) 
n=l  1 


n  =  l 


N 

l  AN“nG  (g) 
n=l  n 


(2.16a) 


(2.16b) 


(2.16c) 


(2 . 16d) 


Substitution  of  Eq  (2.16)  into  Eq  (2.15)  and  factoring  out 
the  constants  aN  and  gives 


V  =  1  +  Vg  +  bNGg 


=  1  +  a..F  +  b  G 

N  N  a  N  a 


(2.17a) 


(2.17b) 


Solving  Eqs  (2.17a)  and  (2.17b)  simultaneously  results  in 


G  -  G  +  A 
_ B _ a _ 

AN(VGa)  +  FBGa 


F  -  F  +  A 


-  A  (F„+G. )  +  F.G.  - 


F  gq 
a  6 


(2.18a) 


(2.18b) 


After  substituting  Eqs  (2.18a)  and  (2.18b)  into 
Eq  (2.8)  and  rearranging  the  terms,  the  final  expression 
in  X  is 


X2N+1 


2N  _  XN+1(  G  }  + 
3  a 


>N(VF'n+1+V<3N  +  1> 


+  X  (F.G  -F  G,) 
6  a  a  6 


^aVW 


+  FN+i(VGe'-Vi(Fa-V  = 


(2.19) 


Equation  (2.19)  is  the  desired  polynomial  expression  in  X 

Note  that  F  ,  F.  ,  G  ,  and  G  are  all  polynomial 
a  6  a  6  c  2 

functions  of  X  .  All  coefficients  of  X  in  Eq  (2.19) 
can  be  calculated,  and  the  2N+1  values  of  X  can  be 
evaluated  numerically. 


Edge  Effects  and  Decentered  Resonators 

When  the  feedback  mirror  of  the  resonator  is 

decentered,  two  equivalent  Fresnel  numbers  can  be  defined 

for  the  cavity.  For  an  aligned  (non-decentered)  resonator, 

Eq  (1.2)  defines  the  N  as 
^  eq 


N  =  ( 

eq 


a2  1  1 

- )  7  (M  -  ±) 

2X  Lg^  2  M 


(2.20) 
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Referring  to  Figure  1-2,  one  equivalent  Fresnel  number  can 
be  defined  for  the  section  of  extending  above  the 

optic  axis,  and  one  for  that  below  the  optic  axis.  Thus, 


N 


eq,U 


[  (1+6) a_] 2 

2X  >Lg]_  2  (M  " 


(2.21) 


=  (1+6) 2  N  „ 

eq 


(2.22) 


Similarly, 


N 


eq,L 


=  ( 1-6 ) 2  N 


eq 


(2.23) 


The  subscripts  U  and  L  refer  to  the  upper  and  lower 
sections  of  . 

Drawing  an  analogy  to  the  aligned  resonator  case,  the 

edge  effects  would  be  expected  to  influence  the  separation 

of  |X.  I  and  I X _  |  when  N  ..  and  N _ _  equal 

1  1'  1  2'  eq,U  eq , L 

3  7 

n  +  g-  or  n  +  ^  ,  n=0,l,2,...  As  6  increases  from  zero 
to  unity,  ^  and  N  L  are  swept  through  a  number  of 

such  points. 

The  values  of  6  for  which  the  edge  effects  should 

alter  the  eigenvalue  separation  are  easily  determined. 

Starting  with  N  T  , 
eq,L 


L 


(2.24) 
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1 


where 


K  = 


0,2,4,...  for  maximum  mode  separation 

1,3,5,...  for  minimum  mode  separation 


(2.25) 


(1-6) 2  N 


eq 


I  +  K 
8  2 


(2.26) 


6  (K)  =  1 


[3  n 

T  +  K 

2  N 

eq  J 


(2.27) 


Equation  (2.27)  is  identical  to  Eqs  (11)  and  (12)  of 
Reference  6. 

Similarly,  for  Neqj[J  , 


N 


eq,U 


|  *  §'+  INT(Neq) 


(2.28) 


where 


0,2,4, 

1,3,5, 


for  maximum  mode  separation 
for  minimum  mode  separation 


(2.29) 


and  where  INT(Ne^)  is  the  integer  portion  of  N 1  .  With  a 

little  rearranging, 


6  (K " )  = 


int (%)*  <i+i 1 


N 


eq 


-  1 


(2.30) 


K  and  K"  can  only  assume  integer  values  in  Eqs  (2.27) 
and  (2.30) . 
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<5  (K)  and  <5(K')  are  specific,  calculated  values  of 


the  variable  6  at  which  the  separation  between  j  A  ^ [  and 

| X _ |  should  be  maximal  and  minimal.  N  T  might  play  the 

dominant  role  in  determining  the  shape  of  the  eigenvalue 

curves  as  functions  of  6  ,  especially  at  large  values  of 

6  .  As  6  approaches  unity,  Neq'L  becomes  much  smaller 

than  N  .  As  noted  earlier,  studies  of  aligned  resona- 

eq,U  y 

tors  show  that  the  depth  of  the  cusps  of  | A  |  decreases  as 
N  is  increased  beyond  N  ,  ...  .  Consequently,  for 

large  decenters,  N  L  may  play  the  dominant  role  in  deter 
mining  the  structure  of  the  eigenvalue  curves. 

The  particular  importance  that  L  has  on  the  shape 

of  the  eigenvalue  curves  has  been  reported  by  Weiner  (Ref  6: 
1830-1).  From  Figure  2  of  Reference  8,  c$  (K)  correlates 
well  with  the  <5  values  at  which  maximum  and  minimum  separa 
tions  occur  between  | A^|  and  | A  |  .  No  correlation 

between  S(k")  and  the  separation  fluctuations  is  made, 
however. 
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III.  Theory  of  Integrated  Intensity  and  Beam  Steering 

A  general  background  theory  of  the  far  field  inte¬ 
grated  intensity  and  beam  steering  is  presented.  General 
relationships  between  the  power  in  a  given  spot  size  and 
the  decenter  parameter  <$  are  discussed.  The  beam  center 
is  defined  and  a  simple  expression  relating  it  to  the  beam 
steering  angle  is  given.  Particular  attention  is  paid  to 
a  Fourier  optics  treatment  of  the  geometric  mode  throughout 
this  chapter. 


Integrated  Intensity 

Integrated  intensity  is  a  measure  of  beam  quality. 

The  integrated  intensity  is  the  amount  of  power  falling 
into  a  given  spot  or  "bucket"  located  symmetrically  about 
the  beam  center.  The  total  beam  power  is  often  normalized 
to  unity.  The  integrated  intensity  or  "power  in  the 
bucket"  is  then  the  percentage  of  total  power  falling  on 
the  spot.  Plots  of  the  far  field  intensity  profile  and  the 
corresponding  integrated  intensity  vs  spot  size  curve  for 
N  =  9.36  and  5  =  0.0  are  given  in  Figures  3-1  and  3-2. 

A  number  of  studies,  both  theoretical  (Refs  6,  12) 
and  experimental  (Refs  13,  14),  have  examined  how  de¬ 
centering  the  feedback  mirror  affects  the  integrated 
intensity.  As  6  is  increased  from  zero  to  unity,  the 
percentage  of  total  power  falling  inside  the  first  Airy 


disk  increases  monotonically  for  the  geometric  mode  (Ref  12) . 
Anen'ev  et  al  (Ref  13)  noted  that  the  axial  brightness  of 
a  Nd: YAG  laser  with  5=1  was  much  higher  than  that  of 
the  same  laser  with  5=0.  The  output  aperture  with 
5=1  was  effectively  twice  as  large  as  that  with  5=0, 
since  the  energy  was  extracted  from  a  single  side  of  the 
feedback  mirror.  Consequently,  the  beam  divergence  angle 
(which  is  proportional  to  the  inverse  of  the  aperture  size) 
was  reduced  and  the  axial  brightness  was  increased. 

The  far  field  intensity  pattern  for  the  geometric 
mode  is  derived  in  Appendix  A.  The  resulting  intensity 
I  (xq)  is 


.  (M-l)  (1  +  5  )  x 

I  (XQ)  =  -rz)  2  t  (M-l)  2  (1  +  5)  2  sine2  [ - - -] 


(M-l)  (1-5  )x 

+  (M-l)  2  (1-5)  2  sine2  [ - pr - 2] 

A  Z 


2  TT  (M+l  )  X 

+  2  (M- 1 )  2  (1-5 2  )  cos  [ - pi - -] 

A  Z 


(M-l)  (1  +  5  )x  (M-l)  (1-5  )x 

•  sine  [ - p> - 2]  sine  [ - P> - 2]}  (3.1) 

A  Z  A  Z 


where  A'  is  the  radiation  wavelength,  z  is  the  separation 
between  the  output  aperture  of  the  resonator  and  the  obser- 


27 


vation  plane,  and  xQ  is  the  coordinate  in  the  observation 
plane.  Setting  M  =  2.0  (corresponding  to  the  cases 
examined  in  this  study)  and  6  =  0.0  ,  Eq  (3.1)  reduces  to 

2  ®7Ixn  x 

r  {xo]  =  n^F[1  +  cos  (irr)]  sin°2  (A)  (3- 

Equation  (3.2)  is  the  far  field  intensity  pattern  for  a 
nondecentered  resonator.  Setting  5  =  1.0  ,  corresponding 
to  a  highly  decentered  resonator,  the  far  field  intensity 
pattern  becomes 

r'  (xo>  =  (xh)2  sinc2  (ri)  (3- 

The  1  +  cos  (6irxQ/ z)  term  in  Eq  (3.2)  causes  the 
energy  to  be  spread  out  more  than  in  Eq  (3.3).  More  power 
will  thus  be  deposited  in  a  given  spot  for  the  decentered 
resonator  (6=1.0)  than  for  the  nondecentered  resonator 
(6=0.0)  .  In  fact,  as  5  is  increased,  the  energy  deposited 
in  the  first  Airy  disk  increases  monotonically  (Fig.  3, 

Ref  6:1832) . 

The  above  discussion  suggests  that  a  highly  decentered 
resonator  may  be  capable  of  depositing  more  energy  into 
a  given  spot  than  a  nondecentered  resonator.  The  actual 
resonator  modes  show  amplitude  and  phase  fluctuations,  while 
the  geometric  mode  has  uniform  intensity  and  phase  profiles. 
However,  since  increasing  6  changes  the  output  aperture 


from  two  slits  to  a  single  wider  slit,  the  resonator  modes 
will  most  likely  show  an  increase  in  integrated  intensity 
for  increasing  6  .  This  is  in  line  with  the  observations 

of  the  previously  cited  works. 

Beam  Steering 

The  beam  steering  is  the  displacement  from  the  optic 
axis  that  the  center  of  the  far  field  pattern  suffers.  A 
knowledge  of  the  beam  steering  properties  of  a  resonator 
is  a  obvious  importance.  Targets  in  the  far  field  might 
be  missed  entirely  if  the  beam  wanders  to  a  great  degree; 
optical  elements  at  the  output  end  of  the  resonator  might 
be  damaged  if  the  beam  is  excessively  displaced  from  the 
optic  axis.  A  particularly  interesting  question  is  how 
decentering  the  feedback  mirror  affects  the  beam  steering. 

The  center  of  the  beam  will  be  defined  as  the  centroid 
of  the  intensity  profile.  For  a  one-dimensional  beam,  half 
the  power  lies  on  either  side  of  the  centroid.  Defining  the 
center  of  the  beam  as  the  centroid  instead  of  the  peak 
intensity  point  can  be  justified  by  considering  Figure  3-3. 
Although  Figure  3-3  is  the  TEM^q  mode  of  a  stable  resonator, 
the  argument  is  the  same.  The  intensity  peaks  are  shifted 
considerably  from  the  cente-  ,z  the  beam,  while  the  centroid 
is  coincident  with  the  centei .  The  centroid  will  in  general 
provide  a  better  approximation  of  the  beam  center;  hence, 
it  shall  be  defined  as  the  center  of  the  beam. 


Fig.  3-3.  Intensity  Profile  of  the  TEM^q  Mode 

Note  that  the  intensity  peaks  are 
displaced  from  the  center  of  the  beam. 

The  beam  steering  angle  0  is  given  by 

x 

0  =  -2. 
z 


where  x  is  the  location  of  the  centroid  in  the  observa 
o 


tion  plane.  Multiplying  by  the  dimensionless  parameter 


2Ma2/A 


2Ma_x 
2  o 


\'z 


2Ma  f 
2  x 


where  f^  is  the  spatial  frequency,  corresponding  to  the 


location  of  the  centroid.  Further,  by  setting 


equal 


to  unity,  the  normalized  beam  steering  angle  becomes 


I 


,  I 
;  I 

/ 

i 


0  .  =  2Mf 

N  x 


(3.6) 


where  the  subscript  N  refers  to  the  normalizations 
employed.  The  half-width  of  the  central  maximum  of  the 
Fraunhofer  pattern  of  a  slit  of  width  2M  is  9  =  1.0  . 

0^  can  thus  be  used  to  relate  the  beam  steering  angles  to 
the  radius  of  the  Airy  disk;  it  will  be  referred  to  as  the 
"normalized  half  angle."  6^  is  a  very  general  parameter, 
as  A'  ,  z  ,  and  a^  are  not  explicitly  required  for 
its  calculation. 

A  major  cause  of  beam  steering  is  tilts  on  the  phase 
of  the  resonator  mode.  If  the  field  immediately  beyond 
the  resonator  output  aperture  is  u(x)  ,  the  corresponding 
far  field  pattern  will  be  k  U(f  )  ,  where  U(f  )  is 

the  Fourier  transform  of  u(x)  and  k  is  a  complex 
function.  A  linear  phase  tilt  may  be  placed  on  u(x)  by 
multiplying  u(x)  by  exp(iwx)  .  Since 

F  [u  (x)  elaX]  =  U(fx  -  w/2tt)  (3.7) 

where  F  is  the  Fourier  transform  operator,  the  centroid 
of  the  far  field  pattern  is  shifted  by  an  amount  propor¬ 
tional  to  the  magnitude  of  the  phase  tilt.  Thus,  tilted 
resonator  modes  will  suffer  beam  steering  proportional  to 
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the  sizes  of  the  tilts. 

Examination  of  Eq  (3.1)  shows  that  the  far  field 
intensity  is  an  even  function  of  xq  ,  regardless  of  M 
or  5  .  The  geometric  mode  will,  therefore,  suffer  no 

beam  steering  at  any  value  of  <5  .  The  phase  (as  well 

as  amplitude)  profiles  of  resonator  modes  show  tilts  and 
aberrations  as  the  resonator  is  decentered.  The  degree  of 
phase  tilting  should  impact  the  relative  beam  steering. 

If  the  resonator  modes  roughly  approximate  the  geometric 
mode,  the  beam  steering  should  be  small.  However,  if  the 
phase  fronts  become  grossly  tilted,  the  beam  steering  will 
be  pronounced. 
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IV.  Error  Analysis  of  FOCAL 


The  computer  code  FOCAL  propagates  the  resonator 
mode  to  the  far  field,  locates  the  bean  centroid,  calcu¬ 
lates  the  integrated  intensity,  and  determines  the  beam 
steering  angle.  The  various  sources  of  numerical  errors 
that  enter  the  calculations  are  examined  in  this  chapter. 
The  errors  are  due  to  the  numerical  techniques  and  approxi¬ 
mations  employed.  As  it  is  difficult  to  quantify  exactly 
the  magnitudes  of  these  errors,  the  following  discussions 
will  refer  only  to  the  relative  sizes  of  the  errors.  A 
listing  of  the  code  FOCAL  may  be  found  in  Appendix  C. 

Far  Field  Intensity  Calculation  Errors 

Calculating  the  far  field  intensity  requires  per¬ 
forming  a  numerical  Fourier  transform  of  the  resonator 
output  field.  Four  primary  sources  of  error  exist  in  the 
calculation.  Errors  can  be  introduced  if  the  mesh  points 
of  the  resonator  mode  are  too  widely  spaced.  Inherent 
inaccuracies  exist  in  the  integration  technique  used  to 
compute  the  Fourier  transforms.  Errors  arise  if  the  far 
field  intensity  is  calculated  at  spatial  frequencies  above 
a  certain  limit.  Finally,  spacing  +-.ie  mesh  points  on  the 
far  field  pattern  too  far  apart  will  create  additional 
discrepancies. 
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If  the  resonator  mode  mesh  spacing  is  too  wide,  high 
spatial  frequency  information  will  be  lost  (see  Figure  4-1) . 

The  field  that  is  calculated  will  be  a  low-pass  filtered, 
and  thus  inaccurate,  representation  of  the  mode.  This  is 
easily  remedied  by  increasing  the  number  of  field  points 
calculated,  with  a  corresponding  decrease  in  the  mesh 
spacing.  Unfortunately,  this  solution  can  become  expensive 
in  computer  time.  A  reasonable  compromise  was  made  in 
this  study  by  calculating  the  field  at  100  points  over  the 
surface  of  the  feedback  mirror.  Increasing  the  number  of 
points  from  100  did  not  produce  a  noticeable  increase  in 
mode  detail.  Consequently,  the  error  introduced  in  the 
calculation  of  the  near  field  pattern  (resonator  mode)  can 
be  considered  negligible. 

Rather  than  use  one  of  the  available  fast  Fourier 
transform  (FFT )  routines,  an  algorithm  based  on  a  numerical 
integration  was  written  to  perform  the  beam  propagation. 

The  integration  technique  used  was  Simpson's  rule  (Ref  15:136-138) 
This  method  connects  adjacent  data  points  with  quadratic 
curves  and  sums  the  resultant  areas.  This  will  produce  some 
error  if  the  data  points  are  not  connected  by  quadratic 
functions  (see  Figure  4-2) .  A  higher-order  scheme,  such  as 
Weddle's  rule,  could  have  been  employed.  However,  this 
would  have  required  additional  computer  time,  so  the  simpler 
(and  less  accurate)  Simpson's  rule  was  used.  This  probably 
introduced  the  greatest  error  in  calculating  the  far  field 
intensity. 


Fig.  4-1.  Plot  Depicting  Improper  Sample  Spacing 
Note  the  loss  of  high  frequency  data. 

The  Nyquist  criterion  states  that  a  signal  with  a 
highest  frequency  component  W  can  be  perfectly  recovered 
if  it  is  sampled  at  a  rate  fg  >  2 VI  and  the  samples  are 
processed  by  a  low-pass  filter  (Ref  16:68-71).  Using  this 
criterion,  the  far  field  intensity  can  only  be  calculated 
for  spatial  frequencies  f  in  the  range 

_  J_  <  f  <  _L_ 

2T  ~  x  "  2T 

s  s 

where  Tg  is  the  mesh  spacing  on  the  resonator  mode.  In 

this  study,  f  was  always  less  than  1/35T  .  As  a 

X  s 

result,  this  source  of  error  can  be  considered  negligible. 


Fig.  4-2.  Exaggerated  Plo,.  Showing  the  Discrepancy 
Between  f  (x)  and  a  Fitted  Quadratic. 

Error  is  Shaded. 

The  final  source  of  error  is  the  spacing  between  data 
points  on  the  calculated  intensity  pattern.  This  is  identi¬ 
cal  to  the  problem  of  mesh  spacing  on  the  resonator  mode  - 
a  wide  spacing  creates  loss  of  higher  frequency  data.  This 
problem  can  be  resolved  by  decreasing  the  data  point  separa¬ 
tion.  During  actual  runs,  the  far  field  intensity  profiles 
were  found  to  be  quite  smooth.  A  mesh  spacing  of  approxi¬ 
mately  0.06  normalized  half  angles  (nha)  was  discovered  to 
be  much  more  than  adequate  to  recover  all  of  the  fine  detail 
The  errors  introduced  in  this  manner  were  thus  quite  small. 

The  use  of  Simpson's  rule  in  the  Fourier  tranform 
routine  introduced  the  most  error  into  the  far  field 
intensity  calculation.  The  other  three  error  sources  were 
small  by  comparison,  and  can  most  probably  be  neglected. 


Errors  Associated  with  Locating  the  Centroid 


The  centroid  of  the  far  field  intensity  is  located 
in  a  three-step  process.  First,  the  total  power  PQ 
between  two  interactively-specified  limits  is  determined 
by  a  numerical  integration.  Then,  starting  at  the  lower 
limit,  the  intensity  is  integrated  until  the  two  mesh 
points  on  either  side  of  the  half-power  point  are  found. 
Finally,  a  linear  interpolation  is  used  to  approximate 
the  location  of  the  centroid  between  the  two  mesh  points. 
This  suggests  three  main  sources  of  error:  part  of  the 
power  is  not  used  in  the  integration  due  to  the  finite 
limits  imposed,  the  exact  power  cannot  be  determined  because 
of  problems  inherent  in  the  numerical  integration,  and  the 
linear  interpolation  will  not  precisely  locate  the  centroid. 

The  far  field  intensity  pattern  theoretically  extends 
from  -w  to  +oo  .  Beyond  a  few  Airy  disks  from  the 
centroid,  though,  the  intensity  becomes  negligible.  Since 
the  numercial  integration  cannot  be  performed  from  -« 
to  +oo  ,  some  of  the  power  is  necessarily  ignored.  This 
introduces  error  into  the  calculation,  as  the  intensity 
profiles  generally  are  not  symmetric.  To  minimize  this 
error,  a  criterion  was  established  that  at  least  90%  of 
the  power  must  lie  between  the  limits  of  integration.  (This 
was  readily  determined.  The  power  in  the  beam  leaving  the 
resonator  was  normalized  to  unity.  As  long  as  the  power 
between  the  limits  of  integration  was  greater  than  0.9,  the 
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90%  criterion  was  fulfilled.)  In  general,  the  power  between 
the  limits  of  integration  ranged  from  90%  to  95%.  This 
was  deemed  to  yield  sufficiently  accurate  results. 

The  use  of  Simpson's  rule  as  the  integration  technique 
introduces  some  error  into  the  calculation.  The  reasons  are 
discussed  earlier  in  this  chapter.  The  far  field  intensity 
patterns  were  quite  smooth,  and  501  mesh  points  were  used 
in  the  calculations  (which  roughly  equates  to  85  points 
across  the  Airy  disk).  Consequently,  the  errors  introduced 
by  the  integration  were  probably  small. 

The  linear  interpolation  adds  further  error  to  the 
calculation.  This  is  illustrated  in  Figure  4-3.  Point  A 
represents  the  actual  centroid.  Point  B  is  the  "centroid" 
located  by  the  interpolation.  An  error  equal  to  e  thus 
exists . 

The  maximum  value  |e|  can  have  is  x  ...  -  x  , 

1  1  n+1  n 

which  is  the  mesh  spacing.  |e|  could  be  reduced  by 
decreasing  the  data  point  spacing.  This  would  require 
additional  mesh  points  and  more  computer  time.  During 
actual  runs,  lx  Ll  -  x  I  was  on  the  order  of  0.025  nha. 
Decreasing  the  spacing  beyond  this  limit  would  have 
required  inordinate  amounts  of  extra  computer  time.  |e| 
could  also  be  decreased  by  using  a  higher-order  inverse 
interpolation  scheme.  Such  a  technique  did  not  exist  as 
a  computer  library  routine,  and  the  time  required  to  create 
a  routine  was  not  deemed  justifiable  in  terms  of  the  poten¬ 
tial  returns. 
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Fig.  4-3.  Linear  Interpolation  Contribution  to 

the  Centroid  Error.  A  is  the  centroid, 
B  is  the  estimated  centroid,  and  e  is 
the  error. 


A  value  of  the  "maximum  error  in  centroid  location 
is  printed  by  FOCAL.  This  is  equal  to  the  mesh  spacing. 

Three  main  sources  of  error  exist  in  locating  the 
centroid.  It  is  somewhat  unclear  as  to  which  introduces 
the  greatest  inaccuracy.  Measures  were  taken  to  reduce 
the  errors  in  all  cases. 

Integrated  Intensity  Calculation  Errors 

The  integrated  intensity  is  calculated  in  a  two-step 


process  after  the  centroid  is  located.  First,  the  intensity 
profile  is  recomputed  between  symmetric  limits  about  the 
centroid.  Then,  using  Simpson's  rule,  the  intensity  is 


integrated  outward  from  the  centroid  to  the  limits.  Errors 
enter  the  calculation  during  the  location  of  the  centroid, 
the  computation  of  the  intensity  pattern,  and  the  outward 
integration.  The  first  two  sources  of  error  have  been 
discussed  earlier. 

The  third  source  of  error  is  again  due  to  the  inherent 
problems  of  the  numerical  integration.  Since  the  intensity 
patterns  were  relatively  smooth  and  the  mesh  spacing  small 
(0.025  nha) ,  the  inaccuracy  of  the  integrations  was  mini¬ 
mized.  Decreasing  the  mesh  spacing  further  would  probably 
have  produced  a  small  increase  in  accuracy  for  the  additional 
expenditure  of  computer  time. 


V .  Study  Results 

This  chapter  presents  the  results  of  the  computer 
analyses  of  the  resonator.  The  effects  of  decenters  on 
the  eigenvalues  are  discussed  first.  Next,  the  changes  in 
beam  quality  due  to  decenters  are  presented.  The  chapter 
is  concluded  by  examining  how  decenters  affect  the  beam 
steering . 

Effects  of  Decenters  on  the  Eigenvalues 

The  eigenvalues  were  studied  as  functions  of  5 

with  N  a  fixed  parameter.  The  twenty-five  cases 

examined  are  listed  in  Table  I.  The  N  values  were 

eq 

chosen  to  lie  about  the  points  of  maximum  and  minimum 
separation  between  |X^|  and  ] X 2 1  (9.37  and  9.87, 

respectively)  and  an  intermediate  value  (9.60).  In  all 
cases,  the  increment  value  for  6  was  0.004.  Typical 
plots  are  shown  in  Figures  5-1  through  5-6. 

A  note  should  be  made  concerning  several  conventions 
used.  6  ( K )  and  <5  00  refer  to  the  <5  values  calcu¬ 
lated  from  Eqs  (2.27)  and  (2.30),  respectively.  6 (K  even) 
refers  to  those  values  of  5 (K)  for  which  K  is  even. 
Similar  remarks  may  be  made  for  6 (K  odd)  ,  6 (K'  even)  , 

and  6 ( K '  odd)  . 

The  eigenvalue  plots  were  analyzed  for  five  different 


items.  First,  the  overall  and  fine  structures  were 


TABLE  I 


Eigenvalue  Analysis  -  Cases  Examined 


Neq 

6  Range 

9.30 

0. 0-0. 4 

9.31 

0. 0-0. 4 

9.32 

0. 0-0.4 

9.33 

0. 0-0.4 

9.34 

0. 0-0. 4 

9.35 

0.0-0. 4 

9.36 

0.0-0. 9 

9.37 

0.0-0. 4 

9.38 

0. 0-0.9 

9.39 

0. 0-0. 4 

9.55 

0.0-0. 4 

9.575 

0.0-0. 4 

N 

eq 

5  Range 

9.60 

0. 0-0.9 

9.625 

0. 0-0.9 

9.65 

0.0-0. 4 

9.80 

0. 0-0.4 

9.81 

0. 0-0.4 

9.82 

0.0-0. 4 

9.83 

0.0-0. 4 

9.84 

0.0-0. 4 

9.85 

0.0-0. 4 

9.86 

0.0-0. 9 

9.87 

0.0-0. 4 

9 .88 

0.0-0. 9 

correlated  to  5 ( K )  and  &(K') 


respectively.  The 


crossing  and  cusping  nature  of  jX^|  and  [ A ^ |  was 
examined  in  detail.  Next,  the  separation  between  |A. 


N  ,  and  <5  . 

eq,L 


and  I A  _  I  was  related  to  N 

1  2 1  eq,U 

The  interleaving  of  the  higher-order  eigenvalues  was 
explored.  Finally,  crossings  between  [X^l  and  [A^ 
were  examined. 

Examination  of  Figures  5-1  through  5-6  reveals 


that  |  A  -^  |  and 


X. 


have  an  overall  periodicity  in 


N  ,  .  This  periodicity  was  observed  in  all  the  cases 
eq,L  r  1 

studied.  The  separation  between  |X^|  and  | X 2 ,  exhibit 

maxima  and  minima  quite  close  to  the  6(K)  values.  In 

the  six  applicable  cases,  6{K=0)  was  somewhat  greater 

than  the  6  value  at  which  the  peak  separation  occurred. 

Superimposed  on  the  overall  oscillatory  structure  is 

a  fine  structure  periodic  in  N  The  structure  is 

eq ,  U 

most  pronounced  for  5  <  0.2  ,  where  N  rT  is  of  the 

eq ,  U 

same  order  as  N 


For  6  >  0.75  ,  the  periodicity 

eq,U  ' 


eq,L 

still  exists,  but  does  not  go  exactly  as  N 

The  separation  between  |X^|  and  [ X 2 !  exhibits 
minor  minima  at  the  <5(K'  odd)  values.  This  agrees 
with  the  theory  of  Chapter  II.  At  these  6  values,  the 
| AjJ  curves  are  depressed  and  the  | \  curves  are 
peaked.  The  magnitudes  of  the  depressions  and  peaks 
decrease  as  6  increases.  The  correlation  between  these 
minor  minima  and  £ ( K '  odd)  is  excellent  for  6  <  0.7  . 
For  6  >  0.7  ,  the  correlation  begins  to  break  down. 
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Fig.  5-1.  Eigenvalue  Plot  for  Ngg  =  9.36,  0.0<6<0.4 
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At  the  <5(K'  even)  values,  two  distinct  types  of 
behavior  are  observed  for  | A ^ ]  and  J  A ^ |  •  In  the 

first  type  of  behavior,  a  minor  separation  peak  is  observed. 
Alternatively,  both  the  |  )  and  1^!  curves  exhibit 

very  steep  slopes  of  opposite  sign.  The  slopes  of  the 
curves  in  these  regions  are  considerably  greater  than  the 
slopes  in  the  surrounding  neighborhoods.  Even  though 
separation  peaks  are  not  exhibited,  the  curves  appear  to 
be  "stretched  apart,"  thus  creating  the  steep  slopes.  This 
form  of  behavior  occurred  more  often  than  the  other.  The 
correlation  between  the  two  types  of  behavior  and  6(K'  odd) 
is  excellent  for  5  <  0.7  .  For  6  >  0.7  ,  the  correlation 
begins  to  break  down. 

Prominent  cusps  or  mode  crossings  between  | A ^ |  and 
j  A 2 I  exist  at  almost  all  6 (K  odd)  values.  Weiner  {Ref  6: 
1831)  asserts  that  at  moderately  large  N  values  (com¬ 
parable  to  those  used  in  this  study) ,  the  first  crossover 
points  do  not  occur  until  5  is  relatively  large.  In 
this  study,  crossings  were  observed  for  5  as  low  as  0.01 

(N  =  9.88,  9.89).  Cusping  can  -  and  usually  does  -  occur 

when  5  is  increased  from  the  first  crossing  point.  Only 
crossing  behavior  is  observed  for  5  >  0.5  . 

The  <5 (k  odd)  cusps  evolve  in  a  regular  manner  as 

Neg  is  increased.  The  four  major  stages  of  the  cyclic 

behavior  are  displayed  in  Figure  5-7.  N  „  and  N  T 

3  eq, U  eq , L 

play  important  roles  in  the  evolutionary  cycle. 
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Fig.  5-7.  Life  Cycle  of  a  Major  Cusp 

(a)  First  Stage?  (b)  Second 
Stage;  (c)  Third  Stage; 

(d)  Fourth  Stage 
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In  the  first  stage  of  the  cycle,  a  deep  cusp  forms 
at  a  pair  of  S(K  odd)  and  5(K'  odd)  values  that  are 
approximately  equal.  Invariably,  6(K  odd)  <  6(K'  odd)  . 

As  N  is  increased,  6(K  odd)  and  &  (K'  odd)  become 
almost  equal,  and  the  cusp  deepens.  In  the  second  stage, 

6(K  odd)  is  . irtually  equal  to  6(K'  odd)  .  The  cusp 
exhibits  two  mode  degeneracy  points.  As  is  further 

increased,  the  third  stage  is  entered.  <5(K  odd)  becomes 
larger  than  5(K'  odd)  ,  and  one  of  the  degenerate  points 
breaks  open.  At  this  stage,  the  cusp  has  evolved  into  a 
single  mode  crossing.  In  the  fourth  stage,  5(K  odd) 
and  6  (K'  odd)  are  somewhat  separated.  The  remaining 
mode  degeneracy  breaks  open,  and  the  eigenvalues  completely 
separate.  As  N  is  increased  further,  the  mode  separa- 
tion  becomes  greater  until  no  cusping  nature  is  evident. 

What  was  once  a  cusp  has  evolved  into  a  region  of  relatively 
large  mode  separation.  This  evolutionary  behavior  was 
observed  in  virtually  all  cusps  corresponding  to  6 (K  odd)<0.4 

The  separation  between  | |  and  | A ^ |  is  a  function 

of  N  TI  ,  N  T  ,  and  6  .  Several  distinctive  types  of 

behavior  may  be  observed. 

If  K  and  K'  are  of  the  same  parity,  several  pairs 

of  6(K)  and  <5(K')  are  approximately  equal,  and  if 

6  <  0.15  ,  then  the  mode  separation  is  characterized  by  rapid 

deep  oscillations.  This  behavior  is  observed  near  N  „  = 
c  eq 

9.37  and  9.87.  In  these  regions,  the  cusps  are  quite  deep 
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and  often  doubly  degenerate.  Mode  crossings  occasionally 
occur.  The  mode  separation  is  quite  unstable  and  changes 
rapidly  for  small  changes  of  5  .  The  behavior  suggests 

that  the  edge  waves  from  both  edges  of  the  feedback  mirror 
reinforce  each  other,  thus  creating  the  large  oscillations 
in  |X.|  and  | X  2 1 

If  K  and  K'  are  of  the  opposite  parity,  several 
pairs  of  6 (K)  and  6(K')  values  are  approximately  equal, 
and  if  6  <  0.15  ,  then  the  mode  separation  remains  roughly 
constant  at  an  intermediate  value.  This  behavior  may  be 
observed  near  N  =  9.625  .  In  these  regions,  changes  in 
6  produce  minor  changes  in  the  separation  between  |X^| 
and  | X 2 |  It  appears  that  the  edge  waves  from  both  edges 

of  the  feedback  mirror  interfere  destructively  with  each 
other.  Each  wave  nulls  out  the  effects  of  the  other. 
Consequently,  the  oscillatory  behavior  is  suppressed,  and 
the  mode  separation  remains  roughly  constant. 

The  greatest  mode  separation  occurs  for  6  =  0.8  , 
near  6 (K  =  0)  .  The  separation  stability  is  quite  good 
near  <5  =  0.8  :  ^  (  |  X ^  |  —  |  A 2  j  )  is  rather  small.  For 
6  >  0.3  ,  the  separation  of  the  peaks  of  the  overall 
oscillatory  structure  increases  with  5  .  A  comparison  of 

the  mode  separations  for  6  =  0.0  and  6  =  0.8  is  given 
in  Table  II.  The  separations  show  improvements  of  38%  to 
51%  when  the  resonator  is  highly  decentered.  This  agrees 
with  the  observations  of  Weiner  (Ref  6:1831). 
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TABLE  II 


Comparison  of  Eigenvalue  Separation 


N 

eq 

6* 

i^l 

|X2I 

AA** 

A%*** 

9.36 

0.0 

1.0515 

0.8107 

0.2408 

42.4 

9.36 

0.776 

1.1437 

0.8007 

0.3430 

9.38 

0.0 

1.0519 

0.8105 

0.2414 

41.7 

9.38 

0.776 

1.1416 

0.7996 

0.3420 

9.60 

0.0 

1.0318 

0.8096 

0.2222 

47.0 

9.60 

0.792 

1.1352 

0.8085 

0.3267 

9.625 

0.0 

1.0272 

0.8079 

0.2193 

51.2 

9.625 

0.788 

1.1363 

0.8048 

0.3315 

9.86 

0.024 

1.0495 

0.8087 

0.2408 

38.4 

9.86 

0.792 

1.1368 

0.8035 

0.3333 

9.88 

0.024 

1.0502 

0.8075 

0.2427 

38.8 

9.88 

0.792 

1.1372 

0.8002 

0.3369 

*  6 

values 

listed  correspond  to  the 

maximum  AA 

values  nearest  6=0.0  or  K=0. 
**  AA  =  IAJ  -  |  A2 1 


***  A% 


6=0 


100 


The  eigenvalue  magnitudes  of  the  higher-order  modes 

interleave  with  and  cross  one  another.  The  interleaving 

is  somewhat  regular  for  6  <  0.6  ,  but  no  periodicity  in 

N ..  ,  N  .  ,  or  6  can  be  readily  discerned.  Some  modes 

eq,u  eq,L  1 

appear  to  meander  with  no  apparent  periodicity  for  6  <  0.6  . 
For  6  >  0.6  ,  the  periodicity  is  much  more  pronounced.  The 
"meandering  modes"  generally  cease  to  exist,  leaving  only  a 
diamond-shaped  interleaving  pattern.  The  pattern  is  roughly 
periodic  in  Ne^  ^  . 

Crossings  between  | X ^ |  and  | X  j  occur  at  all  values 
of  6  .  For  6  <  0.4  ,  the  mode  crossings  are  sporadic  and 
aperiodic.  For  0.4  <  6  <  0.65  ,  the  crossings  are  more 
frequent.  For  6  >  0.65  ,  the  crossings  are  quite  regular 
and  approximately  periodic  in  N  y  . 

In  summary,  the  eigenvalues  exhibit  an  overall  and  a 
fine  structure  periodic  in  N  T  and  N  TT  ,  respectively. 

6 (K)  and  6 (lO  correlate  very  well  to  the  maximum  and 
minimum  separation  points  of  |X^j  and  | X 2 1  .  The  cusps 

at  the  6 (K  odd)  values  show  a  cyclic  behavior  in  N 

eq 

The  separation  between  the  first  two  eigenvalues  is  apparently 
influenced  by  the  edge  waves  from  the  feedback  mirror.  The 
maximum  separation  of  these  two  eigenvalues  occurs  near 
6  =  0.8  ;  the  peaks  in  the  separation  increase  in  magnitude 
as  6  increases.  Finally,  the  higher-order  eigenvalues 
display  a  periodic  interleaving  for  large  6  values. 
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Effects  of  Decentcrs  on  Beam  Quality 


The  far  field  integrated  intensity  of  the  first  four 
modes  was  examined  for  decentered  resonators  with  equivalent 
Fresnel  numbers  of  9.36,  9.625,  and  9.86.  6  ranged  from 

0.0  to  0.9  in  increments  of  0.05.  Specifically,  the  power 
deposited  within  one,  two,  and  three  Airy  disks  (normalized 
half  angles)  of  the  optic  axis  was  calculated.  The  number 
of  Airy  disks  required  to  capture  90%  of  the  power  was  also 
computed.  Plots  of  the  results  for  N  =  9.36  are  shown 
in  Figures  5-8  through  5-15.  Each  plot  shows  the  results 
for  the  geometric  mode  and  either  modes  one  and  two  or  modes 
three  and  four. 

The  plots  were  examined  for  five  specific  items.  First, 
general  increasing  or  decreasing  trends  of  the  data  with 
<5  were  observed.  Second,  the  modes  were  compared  to  determine 
if  any  mode  had  consistently  better  beam  quality  than  the 
others.  The  same  comparisons  were  then  made  between  the  four 
resonator  modes  and  the  geometric  mode.  Next,  the  percent 
differences  of  the  integrated  intensity  for  6  =  0.0  and 
6  =  0.8  were  calculated.  Finally,  the  range  of  6  values 
for  which  the  beam  quality  was  the  best  (most  power  deposited 
in  a  given  spot  size)  was  determined. 

The  integrated  intensity  of  the  geometric  beam  is  dependent 
on  6  .  The  power  deposited  in  one  Airy  disk  increases 
monotonically  as  6  is  increased.  This  agrees  with  results 
published  in  the  literature  (Ref  12).  The  power  in  two  Airy 
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Fig.  5-8.  Integrated  Intensity  in  One  Airy  Disk 

for  Modes  One  and  Two.  N  =9.36 
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Fig.  5-9.  Integrated  Intensity  in  One  Airy  Disk 
for  Modes  Three  and  Four.  N  =  9.36 
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Fig .  5-10.  Integrated  Intensity  in  Two  Airy  Disks 
for  Modes  One  and  Two.  N  =  9.36. 
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Fig.  5-12.  Integrated  Intensity  in  Three  Airy  Disks 
for  Modes  One  and  Two.  Ne^  =  9.36. 
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Fig.  5-14. 


Number  of  Airy  Disks  Required  to  Recover 
90%  of  the  Power  for  Modes  One  and  Two. 


9.36. 
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Fig.  5-15. 


Number  of  Airy  Disks  Required  to  Recover 
90%  of  the  Power  for  Modes  Three  and  Four 


9.36. 
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disks  decreases  slightly  with  increasing  5  .  For  three 

Airy  disks,  the  power  first  increases  until  6  =  0.6  and 
then  decreases.  Finally,  the  number  of  Airy  disks  required 
to  capture  90%  of  the  power  increases  until  6  =  0.2  and 
then  levels  off. 

For  spot  sizes  of  one  and  two  Airy  disks,  all  the  modes 
show  increases  in  integrated  intensity  as  6  increases. 

The  increases  are  more  pronounced  for  the  third  and  fourth 
modes.  The  increasing  trend  is  not  as  obvious  when  the  spot 
size  is  increased  to  three  Airy  disks.  Modes  one  and  two 
display  decreasing  integrated  intensity  for  6  <  0.5  ,  and 
slightly  increasing  integrated  intensity  for  6  >  0.5  .  The 
third  and  fourth  modes  have  slightly  increasing  integrated 
intensity  for  this  spot  size.  The  number  of  Airy  disks  required 
to  capture  90%  of  the  power  generally  decreases  as  6  increases. 
The  first  two  modes  have  widely  scattered  values  as  6  in¬ 
creases,  with  only  a  slight  downward  trend.  Modes  three  and 
four  display  less  scattering  and  have  more  pronounced  decreas¬ 
ing  trends. 

All  modes  have  beam  quality  instabilities  for  0.2  <  6  < 

0.75  .  At  an  instability,  the  beam  quality  deteriorates 
rapidly.  Generally,  the  beam  quality  recovers  from  an  instabil¬ 
ity  within  6  =  ±0.05  .  In  some  instances,  a  series  of 
instabilities  exists,  giving  the  overall  beam  quality  an 
oscillatory  nature.  The  fundamental  mode  always  exhibited  an 
instability  near  6  =  0.7  .  This  instability  has  been  noted 
in  the  literature  (Ref  6:1831-2). 


61 


Individual  instabilities  can  be  traced  through  all 
three  spot  sizes,  and  affect  adversely  the  number  of  Airy 
disks  required  to  capture  90%  of  the  total  power.  For 
example,  the  instability  at  &  =  0.7  for  the  fundamental 
mode  appears  in  Figures  5-8,  5-10,  5-12,  and  5-14.  The 
effect  of  the  instability  on  the  90%  power  point  is  rather 
dramatic. 

No  mode  appears  to  have  better  overall  beam  quality 
than  the  other  modes.  This  can  partially  be  attributed  to 
the  instabilities  that  all  the  modes  suffered.  In  general, 
the  fundamental  mode  has  slightly  worse  beam  quality  than 
the  other  modes.  For  almost  any  specified  6  and  N 
values,  no  method  of  determining  a  priori  which  mode  would 
have  the  highest  integrated  intensity  values  could  be  developed. 

All  modes  have  comparable  or  better  beam  quality  than 
the  geometric  mode  for  spot  sizes  of  one  and  two  Airy  disks 
when  6  >  0.6  .  (This  is  not  true  at  a  beam  quality  instabil¬ 
ity.)  For  a  spot  size  of  three  Airy  disks,  all  modes  have 
somewhat  worse  beam  quality  than  the  geometric  mode,  regard¬ 
less  of  6  .  All  the  modes  require  high  6  values  to  have 

the  same  or  fewer  number  of  Airy  disks  to  reach  the  90% 
power  point  when  compared  to  the  geometric  mode.  Overall, 
at  high  6  values,  the  modes  appear  to  have  comparable  or 
better  beam  quality  than  the  geometric  mode. 

The  percent  difference  A%  in  integrated  intensity  for 
6  =  0.0  and  6  =  0.8  is  defined  as 


A% 


.Power  in  Spot.  _  .Power  j.n  Spot. 

with  6  =  0.8 _ with  5  =  0.0 

'  ~  .Power  in  Spot, 
with  6  =  0.0 


x  100 


(5.1) 


The  percent  differences  were  calculated  for  all  cases  examined 
and  are  listed  in  Table  III.  A  similar  definition  may  be  made 
for  the  percent  difference  in  the  number  of  spots  required  to 
capture  90%  of  the  power: 


A"% 


.Number  of  Spots.  _  Number  of  Spots, 

with  6  =  0.0 _ with  5  -  0.8 

.Number  of  Spots, 
with  6  =  0.0 


(5.2) 


The  A'%  values  are  listed  in  Table  III  under  the  A%  column. 

Table  III  reveals  several  interesting  characteristics. 

In  all  cases,  more  power  is  deposited  in  one  and  two  Airy 
disks  when  6  =  0.8  than  when  5=0.0.  A%  always  monoton- 
ically  decreases  as  the  spot  size  increases.  In  some  instances, 
particularly  for  the  fundamental  mode,  A%  is  less  than  zero 
when  the  spot  size  is  three  Airy  disks.  Modes  two,  three, 
and  four  deposit  90%  of  their  energy  in  smaller  areas  when 
6  =  0.8  than  when  6  =  0.0  .  These  observations  indicate  that 
the  power  densities  near  the  beam  centroids  are  higher  when 
5  =  0.8  than  when  6=0.0. 

In  general,  the  integrated  intensity  valuer  are  the 
highest  for  6  >  0.7  .  The  only  major  beam  quality  instability 
in  this  range  occurs  in  the  fundamental  mode  (•  =  0.7).  This 


l 


63 


(N  CN  CN  <N  CN  CN  <N  (N  <N  CM  CN  CM 


TABLE  III 


Integrated  Intensity  Percent  Differences 


Mode 

N 

eq 

Case* 

A% 

Mode 

N 

eq 

Case* 

A% 

1 

9.  36 

1 

60.5 

3 

9.36 

1 

100. 

1 

9.36 

2 

4.2 

3 

9.36 

2 

22.6 

1 

9.36 

3 

-  7.7 

3 

9.36 

3 

9.5 

1 

9.36 

90 

-55.2 

3 

9.36 

90 

66.9 

1 

9.625 

1 

59.9 

3 

9.625 

1 

93.2 

1 

9.625 

2 

3.6 

3 

9.625 

2 

21.6 

1 

9.625 

3 

-  7.3 

3 

9.625 

3 

7.8 

1 

9.625 

90 

-15.0 

3 

9.625 

90 

69.4 

1 

9.86 

1 

57.9 

3 

9.86 

1 

36.4 

1 

9.86 

2 

6.4 

3 

9.86 

2 

19.7 

1 

9.86 

3 

-  4.6 

3 

9.86 

3 

5.9 

1 

9.86 

90 

6.3 

3 

9.86 

90 

53.1 

2 

9.36 

1 

41.2 

4 

9.36 

1 

95.7 

2 

9.36 

2 

19.9 

4 

9.36 

2 

23.3 

2 

9.36 

3 

9.3 

4 

9.36 

3 

6.6 

2 

9.36 

90 

51.7 

4 

9.36 

90 

33.3 

2 

9.625 

1 

36.1 

4 

9.625 

1 

42.7 

2 

9.625 

2 

13.7 

4 

9.625 

2 

28.8 

2 

9.625 

3 

4.7 

4 

9.625 

3 

7.4 

2 

9.625 

90 

44.9 

4 

9.625 

90 

47.6 

2 

9.86 

1 

81.8 

4 

9.86 

1 

40.5 

2 

9.86 

2 

6.6 

4 

9.86 

2 

17.5 

2 

9.86 

3 

-  2.7 

4 

9.86 

3 

4.0 

2 

9.86 

90 

5.8 

4 

9.86 

90 

27.3 

* 

Case  refers  to  the  number 

of  Airy 

disks  over  which 

the  intensity  was  integrated.  Case  of  90  refers 


to  A"%  values. 
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agrees  with  a  conclusion  from  Chapter  III:  the  beam  quality 
should  be  improved  at  high  values  of  6  . 

Effects  of  Decenters  on  Beam  Steering 

The  far  field  beam  steering  of  the  first  four  modes  was 

examined  for  decentered  resonators  with  N  =9.36,  9.625, 

eq 

and  9.86.  5  ranged  from  0.0  to  0.9  in  increments  of  0.05. 

Additionally,  the  fundamental  mode  was  analyzed  at  a  much 
higher  resolution  for  N  =  9.36  .  For  this  case,  6  varied 
between  0.0  and  0.9  in  increments  of  0.01.  Plots  of  the 
results  for  N  =  9.36  are  presented  in  Figures  5-16  through 
5-18. 

The  beam  steering  data  were  examined  in  three  different 
areas.  First,  the  data  were  analyzed  to  determine  if  the  beam 
steering  angles  had  some  linear  dependence  on  6  .  Next, 

the  oscillatory  nature  of  the  data  was  checked  to  see  if  a 
simple  periodicity  existed,  and  how  the  amplitude  of  the 
oscillations  varied  with  6  .  Finally,  the  means  and 

standard  deviations  of  the  beam  steering  angles  were  computed 
and  compared  to  one  another. 

The  beam  steering  angles  of  the  geometric  mode  are  super¬ 
imposed  on  Figures  5-16  through  5-18.  The  steering  angle  for 
any  6  value  is  either  -0.011  or  -0.01  normalized  half  angles 
(nha) .  Since  the  uncertainty  in  the  location  of  the  centroid 
is  approximately  0.025  nha,  the  optic  axis  lies  within  the 
error  bounds.  This  agrees  with  the  conclusion  reached  in 
Chapter  III  that  the  geometric  mode  should  display  no  beam 
steering,  regardless  of  6 
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Fig.  5-16.  Beam  Steering  Angles  for  Modes 
One  and  Two.  N _  =  9.36. 
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Fig.  5-17.  Beam  Steering  Angles  for  Modes 

Three  and  Four.  N  =9.36. 
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Fig.  5-18.  Beam  Steering  Angles  for  Mode  One 

Under  Higher  Resolution.  N  =9.36 
r  eq 

No  linear  dependence  of  the  beam  steering  angles  on  6 
can  be  discerned  from  the  plots.  A  linear  regression  was 
performed  on  the  data;  the  results  are  listed  in  Table  IV. 
(YlNT '  m'  anc*  r  are  t^ie  Y  axAs  intercept,  slope,  and  correla¬ 
tion  coefficient,  respectively.)  The  magnitudes  of  the 
correlation  coefficients  are  less  than  0.25,  except  in  two 
instances.  This  implies  that  the  data  correlate  poorly  to 
straight  line  fits.  Thus,  the  beam  steering  angles  display 
no  linear  dependence  on  5 

All  the  modes  have  beam  steering  angles  that  oscillate 

pseudorandomly  about  the  optic  axis.  No  period  in  6  , 

N  ,  or  N  T  is  readily  discerned.  The  amplitudes  of 

eq,U  eq , L  1  v 
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TABLE  IV 


Statistical  Analyses  Results  -  Beam  Steering  Angles 


N 

eq 

Mode 

yiNT 

(nha) 

m  r 

(nha/unit  6) 

X 

(nha) 

0 

(nha) 

9.36 

1 

0.0146 

-0.0186 

-0.075 

0.0062 

0.0698 

9.36 

2 

0.0787 

-0.3244 

-0.355 

-0.0673 

0.2565 

9.36 

3 

0.0141 

-0.1403 

-0.205 

-0.049 

0.1924 

9.36 

4 

-0.1935 

0.2294 

0.133 

-0.0902 

0.4865 

9.625 

1 

-0.0017 

-0.0061 

-0.03 

-0.0044 

0.0576 

9.625 

2 

0.0256 

-0.225 

-0.24 

-0.0757 

0.2632 

9.625 

3 

0.0677 

-0.1729 

-0.271 

-0.0101 

0.1794 

9.625 

4 

-0.0106 

0.0943 

0.095 

0.0318 

0.2781 

9.86 

1 

-0.0162 

0.0075 

0.041 

-0.0128 

0.0516 

9.86 

2 

-0.0261 

-0.0507 

-0.09 

-0.0489 

0.1583 

9.86 

3 

-0.1311 

0.0585 

0.055 

-0.1048 

0.3016 

9.86 

4 

-0.0181 

-0.0327 

-0.033 

-0.0328 

0.28 

the  oscillations  appear  to  be  random.  The  oscillatory  nature 

is  more  apparent  in  the  high  resolution  case  (Figure  5-18) . 

Again,  the  period  does  not  correlate  to  8  ,  NeQ  ^  ,  or 

N  .  .  It  is  interesting  to  note  that  the  beam  steering 
eq,L 

instability  at  6  =■  0.71  is  coincident  with  a  major  beam 
quality  instability  (see  Figure  5-8). 
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The  mean  values  x  and  the  standard  deviations  a  of 
the  beam  steering  angles  were  computed  for  all  cases.  The 
results  are  listed  in  Table  IV.  The  mean  value  is  always 
less  than  0.105  nha;  in  two-thirds  of  the  cases,  it  is  less 
than  0.05  nha.  A  rather  simplistic  error  bound  on  the  mean 
values  is  0.025  nha.  Thus,  the  mean  beam  steering  angles  are 
quite  close  to  the  optic  axis.  The  standard  deviations  are 
generally  less  than  0.3  nha.  No  standard  deviation  exceeds 
0.5  nha.  The  fundamental  mode  quite  clearly  has  the  lowest 
standard  deviation.  Its  standard  deviation  is  approximately 
one-third  (or  less)  of  the  standard  deviations  of  the  other 

modes  for  any  N  „  value. 

1  eq 

In  summary,  the  beam  steering  angles  oscillate  somewhat 
randomly  about  the  optic  axis.  No  linear  dependence  of  the 
beam  steering  angles  on  6  is  apparent.  The  average  values 
of  the  beam  steering  angles  are  quite  close  to  zero,  and  the 
standard  deviations  are  generally  less  than  0.3  nha.  Finally, 
the  fundamental  mode  has  the  lowest  beam  steering;  its  standard 
deviation  is  one-third  of  or  less  than  the  standard  deviations 
of  the  higher-order  modes. 


VI .  Unstable  Resonator  Design  Criteria 

This  chapter  presents  several  simplified  design  criteria, 
based  on  the  observations  described  in  Chapter  V.  Design 
criteria  for  decentered  resonators  are  presented  first,  followed 
by  criteria  for  nondecentered  resonators.  It  is  important  to 
realize  that  the  comments  below  are  based  on  properties  of 
resonators  with  9.3  <  N  <  9.9  .  The  criteria  may  vary  for 
resonators  with  N  values  removed  from  this  range. 

Decentered  Resonator  Design  Criteria 

From  the  observations  of  Chapter  V,  highly  decentered 
resonators  offer  certain  advantages  over  nondecentered 
resonators.  The  beam  quality  tends  to  be  better  and  the  mode 
separation  greater  for  highly  decentered  resonators.  A 
cavity  can  be  designed  to  exploit  these  advantages.  The 
following  guidelines  are  advanced  for  decentered  resonator 
designs . 

(1)  The  resonator  should  be  designed  for  operation  in 
the  fundamental  mode.  This  follows  directly  from 
the  beam  steering  results.  The  fundamental  mode 
has  considerably  lower  beam  steering  angles  than 
the  higher-order  modes.  The  standard  deviations 
of  the  fundamental  mode  beam  steering  are  approxi¬ 
mately  one-third  c  or  less  than  the  standard 
deviations  of  the  higher-order  modes.  The  mean 


70 


beam  steering  angles  for  the  fundamental  mode 
are  extremely  close  to  the  optic  axis.  Operation 
in  the  fundamental  mode  will  generally  yield  the 
lowest  beam  steering  angles. 

(2)  The  resonator  should  be  designed  for  very  high 
decenters  (near  5  (K=0)  )  .  For  the  particular  N 
range  studied,  5  should  be  approximately  0.8. 

This  value  is  based  on  the  mode  separation  and  the 
separation  stability,  the  integrated  intensity,  and 
mechanical  design  considerations. 

The  greatest  separation  between  |A.|  and  the  magnitudes 
of  the  higher-order  eigenvalues  occurs  for  6=0.8  .  This 
value  of  6  will  enhance  the  single  mode  operation  of  the 
resonator.  The  separation  is  quite  stable  for  5  =  0.8  ;  only 
small  changes  in  the  separation  will  occur  if  6  deviates 
slightly  from  the  operating  point. 

Highly  decentered  resonators  tend  to  deposit  more  energy 
in  a  given  spot  size  than  nondecentered  resonators.  The 
percent  difference  values  recorded  in  Table  III  illustrate 
this  point.  The  fundamental  mode  deposits  more  energy  in  spot 
sizes  of  one  and  two  Airy  disks  when  5  =  0.8  than  when 
<5  =  0.0  .  The  overall  results  indicate  that  the  power  densities 
near  the  beam  centroid  are  greater  when  5  =  0.8  than  when 
<5  =  0.0  for  the  fundamental  mode.  The  higher-order  modes 
also  show  beam  quality  improvements  at  high  decenters.  Any 
higher-order  modes  existing  in  the  resonator  would  thus  have 
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improved  beam  quality  when  6  =  0.8  . 


The  mechanical  construction  is  simplified  if  the  reson¬ 
ator  is  highly  decentered.  Nondecentered  resonators  typically 
must  be  built  with  the  feedback  mirror  mounted  on  a  transmissive 
optical  element  or  suspended  by  a  spyder.  Consequently,  the 
outcoupled  beam  is  partially  blocked  by  a  mechanical  element. 
This  can  be  particularly  bothersome  if  the  power  densities  are 
very  high.  (Krupke  and  Sooy  (Ref  17:580-1)  have  reported  the 
use  of  an  annular  scraper  mirror  to  avoid  these  difficulties.) 

In  a  highly  decentered  resonator,  the  spyder  sizes  can  be 
reduced  and  the  clear  aperture  size  increased.  Notably,  if 
6  =  1.0  ,  no  spyder  or  other  interfering  mechanical  mounting 
device  need  block  any  portion  of  the  outcoupled  beam.  This 
is  illustrated  in  Figure  6-1. 

Nondecentered  Resonator  Design  Criteria 

The  results  of  the  eigenvalue  study  suggest  an  optimum 
equivalent  Fresnel  number  for  the  resonator  if  fundamental 
mode  operation  is  desired.  This  value  is 


N 


eq , opt 


n  5/8 


n=0 , 1 , 2 ,  .  .  . 


(6.1 


N  .  results  in  the  best  mode  separation  stability.  Small 

accidental  misalignments  of  the  resonator  mirrors  will  have 
little  effect  on  the  mode  separation  at  this  value  of  N 

eq 

Essentially,  the  edge  waves  from  the  feedback  mirror  tend  to 


OUTPUT  APERTURES 


NEAR  FIELD  BEAMS 


1 


null  each  other  out.  The  result  is  a  relatively  constant 
mode  separation  for  6  <  0.15  .  Perkins  and  Cason  (Ref  18: 

200)  report  a  similar  choice  for  N  ,  ,  but  do  not  justify 

its  choice  through  edge  effects  arguments. 

Limitations 

As  noted  in  the  introduction ,  the  above  design  criteria 

are  based  on  studies  of  resonators  with  9.3  <  N  <  9.9  . 

eq 

The  criteria  may  not  be  particularly  applicable  for  resonators 
with  very  high  or  low  equivalent  Fresnel  numbers.  Any  specific 
design  may  deviate  from  these  criteria  because  of  effects 
not  analyzed  or  considerations  not  taken  into  account  in  the 
above  discussions. 
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VII.  Conclusions  and  Recommendations 


The  conclusions  of  this  study  are  summarized  in  this 
chapter.  Areas  that  warrant  additional  research  are 
outlined . 


Conclusions 


(1)  Two  equivalent  Fresnel  numbers,  N  ..  and  N  T  , 

eq,U  eq,L 

can  be  defined  for  aecentered  resonators.  N  ,,  and  N 

eq,U  eq,.L, 

are  functions  of  the  feedback  mirror  fractional  decenter  5  . 
The  defining  relationships  for  the  two  Fresnel  numbers  are 
given  in  Eqs  (2.22)  and  (2.23). 

(2)  Drawing  an  analogy  to  the  case  of  a  nondecentered 
resonator,  two  functions,  5(K)  and  6  (lO,  can  be  derived 
that  yield  5  values  at  which  the  separation  between  jx^| 
and  | A 2 |  is  minimal  or  maximal.  |X  |  and  ] A 2 |  are  the 
magnitudes  of  the  two  lowest-order  eigenvalues.  K  and  K " 
are  integer  variables  used  in  the  functions  defining  6 ( K ) 
and  6 (K ' )  ,  respectively  (see  Eqs  (2.27)  and  (2.30)). 

(3)  The  eigenvalues  are  functions  of  N  „  and  N  T  . 

eq , U  eq,L 

Two  types  of  structure  are  exhibited  in  the  } A ^ [  and  | A 2 I 
curves:  an  overall  structure  periodic  in  N  T  ,  and  a 

superimposed  fine  structure  periodic  in  N  ^  .  6(K)  and 

6(K')  quite  accurately  predict  the  values  of  5  at  which  the 
separation  between  | A ^ |  and  | A 2 i  is  maximal  or  minimal. 

(4)  Major  cusps  in  the  |A^|  and  | A 2 1  curves  occur 
at  the  5 ( K  odd)  values.  The  cusps  undergo  a  regular 
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is  increased: 


evolutionary  pattern  with  four  phases  as  N 

a  very  deep  cusp,  a  doubly  degenerate  cusp,  a  mode  crossing  with 

a  single  degenerate  point,  and  a  complete  separation  of  the 

modes.  The  cyclic  behavior  is  apparently  related  to  N  ..  , 

N  T  ,  and  constructive  and  destructive  edge  effects  from 
eq,L 

the  feedback  mirror. 

(5)  The  mode  separation  maxima  between  |A.  |  and  1^1 
at  the  6 (K  even)  values  tend  to  increase  in  magnitude  as 

6  is  increased.  The  greatest  separation  invariably  occurs  when 
5  =  6 (K=0)  .  At  this  6  value,  the  mode  separation  is  very 

stable . 

(6)  The  higher-order  eigenvalues  interleave  somewhat  ran¬ 
domly  at  low  6  values.  The  interleaving  becomes  regular  and 
roughly  periodic  in  N  ^  for  6  >  0.65  .  The  crossings 
between  j  A ^  I  an<3  U3|  behave  similarly. 

(7)  The  power  deposited  in  one  and  two  Airy  disks  gener¬ 
ally  increases  as  5  increases  for  all  modes.  The  trend  is 
somewhat  nebulous  when  the  spot  size  is  increased  to  three 
Airy  disks.  The  number  of  Airy  disks  required  to  receive  90% 
of  the  power  generally  decreases  as  6  increases.  This  trend 
is  somewhat  vague  for  the  fundamental  mode,  but  is  quite  pro¬ 
nounced  for  the  third  and  fourth  modes. 

(8)  All  modes  display  beam  quality  instabilities  for 
0.2  <  6  <  0.75  . 

(9)  No  mode  consistently  deposits  more  energy  into  a  given 
spot  size  than  the  other  modes.  On  the  average,  the  fundamental 


mode  appears  to  have  slightly  lower  beam  quality  than  the 
other  modes. 


(10)  The  beam  steering  angles  have  no  linear  dependence 

on  6  .  This  is  confirmed  by  the  low  correlation  coefficients 

from  the  linear  regression  analysis. 

(11)  The  beam  steering  angles  for  all  modes  fluctuate 
pseudorandomly  about  the  optic  axis.  No  periodicity  in  6  , 

N  U  '  °r  Neq  L  reac^ily  apparent.  The  amplitudes  of 

the  fluctuations  display  no  general  trends  with  6  ,  N 


or  N 


eq ,  U  ' 

eq,L 

(12)  The  means  and  standard  deviations  of  the  beam  steer¬ 


ing  angles  are  lowest  for  the  fundamental  mode.  The  standard 
deviations  of  this  mode  are  at  least  one-third  as  small  as  those 
of  the  higher-order  modes.  The  standard  deviations  of  all  modes 
tend  to  be  less  than  0.3  normalized  half  angles,  although  in 
one  case,  6  =  0.5  normalized  half  angles. 


Recommendations 

(1)  The  analyses  of  this  study  were  performed  on  a  limited 
class  of  resonators  (9.3  <  Necf  <  9.9,  cavity  magnification  =  2). 
The  same  analyses  should  be  performed  on  resonators  of  lower 

and  higher  equivalent  Fresnel  numbers  and  various  magnifications. 
It  should  be  determined  if  the  observations  and  results  of  this 
study  apply  to  other  classes  of  unstable  resonators. 

(2)  The  beam  steering  angles  should  be  correlated  to  the 
phase  tilts  on  the  output  modes.  The  tilts  might  have  to  be 
weighted  to  account  for  intensity  fluctuations  and  the  asymmetric 
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shape  of  the  output  aperture. 


(3)  The  beam  quality  instabilities  might  correlate 

to  large  resonator  mode  aberrations.  This  possibility  should 
be  examined,  and  a  model  developed  to  predict  at  what  6 
values  the  instabilities  will  occur. 

(4)  As  noted  in  Chapter  V,  a  beam  quality  instability 

and  a  large  beam  steering  angle  exist  in  the  fundamental  mode 

for  6  =  0.71  ,  N  =  9.36  .  This  observation  should  be 

eq 

investigated  further  to  determine  if  there  is  a  general  cor¬ 
relation  between  beam  quality  instabilities  and  beam  steering, 
or  if  this  is  merely  an  isolated  coincidence. 
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APPENDIX  A 


Far  Field  Intensity  of  the  Geometric  Mode 

The  far  field  intensity  pattern  for  the  geometric  mode 
of  a  decentered  resonator  is  derived  in  detail  below.  The 
geometric  mode  has  a  uniform  amplitude  and  phase  on  any  plane 
perpendicular  to  the  optic  axis  inside  the  resonator  (see 
Chapter  I).  The  mode  amplitude  is  exactly  zero  beyond  the 
shadow  boundaries.  In  this  derivation,  the  amplitude  is 
arbitrarily  set  equal  to  unity  and  the  phase  to  zero. 

Referring  to  Figure  A-l,  the  transmission  function  of 
the  output  aperture  of  the  resonator  is 


t  (x ) 


=  rect[ 


x  -  >5  (M+ 1 )  (1  +  6)  . 
(M-l)  (1  +  6)  1 


+ 


rect[ 


x  +  h  (M+l)  (1-6) , 
(M-l)  (1-5)  J 


(A .  1 ) 


where  a2  has  been  normalized  to  unity  and 

r 


rect  (x) 


1  M<J5 

0  elsewhere 


Assuming  that 


z 


>  > 


kM2 
2  ' 


(A. 2) 


(A. 3) 


the  Fraunhofer  diffraction  formula  may  be  used.  Here,  z  is 
the  distance  to  the  far  field  observation  plane  and 
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k  =  2 tt/ A'  is  the  propagation  constant.  Then  (Ref  8:61), 


ikxo 

exp  (ikz)  exp  (--■  ■) 

U(xQ)  =  - - — —  F[  U  (x)  t  (x)  ] 


(A. 4) 


f  =  T-T— - 
X  A  Z 


where  F  is  the  Fourier  transform  operator,  u(xq)  is  the 
far  field  distribution,  and  U(x)  is  the  geometric  mode. 
Performing  the  indicated  operations  yields 

ikXQ 

exp  (ikz  )  exp  (-  -5--°) 

n  \  _  Z  Z 


•  {  (M-l)  (1  +  5)  exp[ -imf  (M+l)  (1  +  5)]  sine  y 

X 


+  (M-l)  (l-6)exp[  iTtfx(M+l)  (1-5)]  sine  n)  (A. 5) 


where 

y  =  (M-l)  d  +  <5)fx  (A. 6a) 

n  =  (M-l)  (1  —  6 )  f  x  (A. 6b) 

The  far  field  intensity  is  of  primary  importance;  it  is 

given  by  the  product  of  U(xQ)  and  its  complex  conjugate 

* 

U  (xq)  .  Denoting  the  intensity  as  I  (x^)  ,  after  some 
straightforward  manipulations  I  (xq)  reduces  to 
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Fig.  A-l .  Output  Aperture  of  the  Resonator 
Note  a2  has  been  normalized  to  unity. 


APPENDIX  B 


Program  Listing  -  EIGEN 


The  computer  code  EIGEN  calculates  and  plots  the  eigen¬ 
values  of  a  decenter ed,  unstable  resonator.  Two  modes  of 

operation  are  allowed:  N  can  be  held  fixed  and  6  varied, 

eg 

or  5  can  be  held  fixed  and  N  varied.  Tables  of  the 

eq 

eigenvalues  are  generated  and  written  to  TAPE?.  The  program 
is  essentially  a  modified  version  of  the  eigenvalue  routine  of 
BARC2  (Ref  3:61-74).  The  modifications  include  the  provisions 
for  iterated  calculations  and  the  addition  of  a  plotting 
routine . 

The  required  inputs  are  listed  below.  All  inputs  are 

real  variables,  except  as  noted. 

IRES:  flag  -  enter  N  if  N  is  the  variable 

eq 

and  any  other  letter  if  i  is  the  variable 
vi\G:  cavity  magnification 

The  following  inputs  are  required  if  N  is  the 
variable : 


NEQMIN:  minimum  N  value 

eq 

NEQMAX :  maximum  N  value 

eq 

NUMNEQ:  number  of  N  values  at  which  the 

eq 

eigenvalues  are  to  be  calculated 
DELMIN:  6  value  (fixed  parameter) 

The  following  inputs  are  required  if  6  is  the 
variable : 
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DELMIN:  minimum  6  value 

DELM^X :  maximum  6  value 

NUMDEL:  number  of  <S  values  at  which  the  eigen¬ 

values  are  to  be  calculated 

NEQMIN:  N  value  (fixed  parameter) 

eg 

The  IMSL  mathematics  library  and  DISSPLA  plotting 
package  are  required  for  program  execution.  Since  DISSPLA 
is  used,  the  code  should  be  submitted  as  a  batch  job. 
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IF<  l  AE  J.NE.INEQIC AL  L  G  AF  <  i)t  LP  I  *.  *  'jH  -C  AL  t  «  Of  L  *  AX  »  .  i>  *  >  *1  .  C  AL  L  *  1  .  1  *  I 
OUTPUT  AUN  DATA 

IFI  :Rr-,#C0.!f.Ea>M*  I  ^  <  7  ,2u' .  >M  'JPl\*LEUF  AX.fvEQT  \C  •  Q  *0!  t  PI  * 

IFC  l*es.NE.  If.tu  >Ur  ;  :t  <  7*201''  >Otl  «r.  •ottf  AX.OfU  -.C#  vuHOf  L»M  o-l  \ 
2000  F0AMAT<1M1.1X.*£  IGEmWAlUl  TAHll  AI  A  FtNCTt**.  F  *.t  Q*  //I X,  •PA\f.f 
iof  *ta:  •*F*.4#ix*«  ti  •*  :x*f  -  .4/lx  *•  .i  i  :ncfim:*t  valik  =  •  *rs.i>/i» 

2***U»1Hf<  OF  fiFQ  WAL'Jr  j  =  •  *  t  2/1  •  Ot.L  T  A  CC0V3  T  A  •  T  »  -  **F'.‘» 

2010  FOA*AT<  1H1.1X  *•£  ICt*  LVAlU*  fAULt  >  A.  A  FLNCTITN  .'F  D:LTA*//l***n  ANG 
It  Jf  DELTA:  •  «F>»t.t  |X«*~U*«lXtt  .r/1 *, • U£ LT A  !  \  C  T :  *  t  N  T  VALUE  =  *.F 
2a.6/lX**NUMHEP  OF  OLlTA  V  AL  UL  b  -  •  *  l  3/ 1  X*  *  *.t  0  lC0Sr,TANr>  =  •#F*?.  4) 

MSUPN<I)  =  MAG*  •  <  I  -  1  > 

HSUBHII)  =  1  ♦  1/MAG**2  ♦  ♦  1  /HAG*  •  <2*1  -2 1 

HSUtMIl  1  =  1. C 

MSUPNC  i >  =  i 

DC  10  [=2*51 

MJUPM  I  »  =  MAG*MSUPM  1-1  ) 

MSUUNC  )  =  NoUUN(I-l  >♦!  ./MLUP.iCl  »**2 
10  CONTINUE 

RNfUG  =  AL0G<250.*Nf  UMAX  )  /  ALwGi  MAG  I 

I  F<  P.AiHl  G.GC  .5C  •  )GL  M  7  77 

NUI G=  INF <PSHIG) *1 

MAI  T  E ( 7  *9  9b  )RNB:G#  «BIG 

996  F  0**  MA  f  <  1  <#•  CALCUL  A  HO  N.jIG  =  •  *fi  /IX  ,*Sf  LECT£  Q  SHIG  -  •*!?///> 
FIMSH  1M1  T  £  ALI  2AT  I  LNS  A%0  jET  UP  00  L-"»P'- 
NEO  =  PiF  1  MI  N 

00  20  2  0  1  1MQ=1  tf.UMf.l  U 
T  =2 • *NC  Q*PI *MAG*MAG/ < FAo*PAG-i  . > 

DEL  T  A=OLLM:  f* 

00  2 030  IOELTA=l tNUMOEL 


t 
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COMPUTE  COrfFICtfM  j  OF  T  Ht  roL  YNl  AL 
Pil\  =  COtf <1  I* /«  *  .Jt  G  ♦  C-l>  iJ  <NOtG-i  »  ♦  ... 

♦  COCFtNOEGW  ♦  CO LF  !M)LG*1  I 

ALPHA:-!. ♦OFLT A 
BCT  A  = l* ♦DEL T  A 
COEFlil-CHPLXU.fj.  I 
HOlG  =  2*.BIG*1 
NCJ£F=f.O£G»l 
N^m6«l 
O'.'  15  l=i*“ 

A  'U  =CSOiUI<i.*£Yt:*PI  •  T/f’jUBMZM 
A ti2  =-T*k:Y£/HSUH\!  I  ) 

A  S3  =Bl  r  a*  u  .-1  ./h  , u»  M  :  *1  ) ) 

AM  =  ALP  HA  •  !  l.-l./HiUP'.OlU 
AN6 -Ht  f  A-ALFHA/'-  SGPM  I  M  I 

ano  =  alpha-b»  ta/pdop*.!  i 

FOt  TACI  I  =-! CCXP  !  Af\2  •  AN  3 •• 2  I / AN  51/A  » 1 
F ALPHA!  II=*!C£XF!  ANL-AV.**.  I/AV  I /Ai.l 
GHFTA! i  IMCk.XP!  Af,2*  A  Lb  -  *0  )  /  A  r.fe  *  /A,.  I 
GALPHAC!  I  =!C£Xf <AN2  *AM  «*2 1/  AM  i/AM 
15  CONTINUE 

F  ALPHA!  HI  =CFXP!  AN2«  Ml  T  A  •  •  2 I  /  Hi  T  A/ (  -  AM  I 
Fb£  IA(MI=FAlPHA(h  ) 

GHf  T  A(M)=CI  XP!  A\2*  ALPHA  •  •2>/ALr'H  A /AM 
GAL  PH  A!  M>  =GBL7 A(*> 

C C-L F  !  2  I  =-  !f  b£T  A  !  1  )♦  b ALPHA  ! 1  1*1,1 
Ll  =  *.  BIG*2 
%  A:  1 
NB=  l 

DC  21  1=3*11 
X  2=  CHPL  X! 0. ,0.1 

Y  2=  X  2 

DC  Id  JA= 1  *  *  A 
KA=  A-JAM 

XUFHLTA!  JA  )*GALPhA  (K  A  I  -F  ALf  hA  (  J A I •  Gb£  t  A !  KA  I 
X2  =  XI  ♦  X  2 

18  CONTINUE 
IFU.ea.3l  GO  TC  20 
O';  1*!  JH=l*f.B 
KH=’.H-JU»  l 

Y  1  =  F  ALPHA!  Obi  •GHL  TA  (Kbl-FBL  :  A!  Jtl  I  •  GALPHA  (  KB  I 
Y2=Y1 *Y2 

19  continue 

Nil:  uH  ♦  1 

20  2I=FBtf A! 1-2 1-FDtt  AU-i  I 

2 ?  =  CALPHA  < ! -2 l-CAL  ° mA ( I -1  I 
C Oc F 1 1  >  =  K2*Y2'*Z 1  ♦Z.’ 

N  As  N  A  ♦  1 

21  CO*?IN»JF. 

L2=.BI G*5 
N A=  BIG-1 
NH=.;BIG-1 

DO  2d  I  =L2,  f.CCFF 
X J  =  CMPLX! 0. , 0 • ) 

Y  2  =  X2 

if! i.eo.scaFi  bu  r ?  23 

00  2 2  JA  = 1  *  f  A 
XA=VrtrG*JA-/.A 

X  1  =  FHk  T  A(H-JA)»GALP^A!KA  l-FALPHA!  *-JAI  •GHl  T  A  !  K  A  I 
X  2=  Xl  ♦  X2 

22  c^rr.ut 

IF!  ;.eo.L2i  oc  7o  2  •- 

23  OC  :•*  JH=  t ,  r  H 
K  B=  .U [ G* JH- '  H 

Y 1  =  F  AL  pHA  ! K- JH I  •  GHL  r  A  (Kill  *F  Mt  T  A  (  *  - I  »GAl  PH 4  !  K  t) 


#* 
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T2=u*r? 

24  CwNTINUE 
NH=  .tt-i 
60  IJ  21 

25  00  2a  JH- l •  LB*  G 

Kt»= 

r J =F  AL^MAi JHI -GBE T A IKdl -f  HE  fA(JH>*GAL^hA(KH) 

¥ ?=  ¥!♦¥.? 

2fc  CONTINUE 

27  /IsFHt  TACMI«<CALPMA  <  I-*- 1  |-GiU  :  AC-M-l  1 1 

22=  teALPHA  <**)•<*  AlPHAC  f  All-H-i  J> 

ecu  <  ii-x^yjww? 

N  A=  '*  A-  1 
CON  Tif.UL 

COMPUTE  HJCTS  Lf  fJLYN>IAL  ^Hd  IH'.L  ZCP'LV*  T*. 

OBTAIN  rrtC  f  IGLNVAl  Ut  j  t  AND  Thl  \  w-Otr  1 1GE  NVALLf  S  BY 

size. 

CALL  /CPOLYKGU  tM)‘  G#LAMMUA  »It  *1 
UK  I  TC(7#o  nMQiDUr  A 
89  FGalMAT(lX»liC(lH*)//lX«  *M*G  "  •  «F  -»  .4 /l  x  »  •  OL  L  T  A  =  •  fFK.6//lX* 
l?  X*  *L  AMtJOAC  KtAL  !•  #3  *♦  *L  AMHO  A  (  i  ** Au  I  •  •  o  X  *  •  £ VMAC*  ,  1  3  X  ,  *C  VPri*/  ) 

I  =1 

DO  7  J  I1=2*».0EG 

SI/C'  KCALCLAMbOAU  II  •  •  C  ♦  A I  m  A  o  <  L  A  MtiU  A  C  *  I  I  *  ♦  2 
K  =1 

do  ib  j=u#f.oec 

SI/ El  =-u:AL<  LAMBDA  (J  M  *  •  -♦  AI  M  A  G<  l  A  Kb  0  A  <  J)  1**2 
IF(  SIZE  l.LT  .SIZE  >  Gti  TC  75 
K  =  J 

SIZCsCI/El 
75  CON  1 1  HUE 

C0JH  =  L  AMilOACI  I 
LAHHOA(  l  >=LAMHOA(M 
LAMBOA(K) =COUM 
CL(  I)=LAHrtDA<  II 

EVPM5A1  A..2C  A  I  HAG  (CL  C I  I  >  ,KfAL«CLC  III-IO./PI 
SMA=?EAL<CL(l I  I  •  AIMAutCc  C  Z  >  I  *«2 
SHAO  =  SQKT<  MAI 

URI  T£( 7*333)1  .LAHHOAd  I.SMAGtlVFH 
333  F0R.MATUX«IS»lXt4<G:4.7tlXM 
I  =1  1 

70  CON  TIN Jt 

CVpM=AfANJ<AlMAG<LAvOi)A(NiU  Gl  I  •  • E Al ( LA MdOA < LJEG I  I  ) • la J . /PI 
SmA  =9EAl( LAMBDA (NOE L> I  *  *2 #A I MAu< l AMU D AC  SDEG  1 1  * • ? 

SMAG  =  $GKT  (  MAI 

Uhl  TL  (  7*3J3INDCGfLAKHOA(r«OLGI*.MAG*EVPH 
C 

C  SEARCH  FOR  MAX  AND  M I N  MODF  SEPARATION-*  AND  MAX  AND  MIN 

C  EIGENVALUES. 

C 

I  FT  1 1  NED#  I OFLTA-3  USAu  #2050  *  20*>. 

2040  MAG12=CAJ0(LAMUDA( l  II 
MAG?2  =  CAB »f LAMBDA (2  1 1 
SfP?=MAG12-MAC22 
60  TO  20  70 

2050  LAMIU  ULAMHDACt I 
LAM 021  =lAMHOAC  I 
LOC  I  =CHPl  X (  M  3,  Jtl  A) 

MAC  1 1  =C  A*J  >T  LAMBDA  U  II 
MAG  ?  I  =  C  Art  ,U  AMHO  AC.'  1 1 
3CP1 =MAC1 l-MAGCl 
CD  TO  2C7t 

2060  MACIJ=CAH j(l AMH0A< 1 II 
MAG  2  ) =C  ArtST l  A**tiOA  (  2  >1 
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SEPg-MAG10-MAG2u 

Iff  MAGI 1  •  GT •MAG12«A  <0«MAGll • GT •  MAGI L ) U  » i IL < 7 ,2 ?  • 3  U AMttl  1 #MAGll*LrC 

II 

IFCMAC1  l.U  .MAG12.A  O.MAiiU  .L7.“AGK  >y»  i TEC  7*2; -3  TLAMHll  **AGll«L'C 

It 

IFCMAG21.67  .HAG^.A’.O.^AGri  .GT,kaG2j>*'  JTfl7*2122>.  A*H21  »  MAG21  •  L  C 
11 

lfi  MAG21.LT  .MAG22.A  .0. “ AG  21 . L T . MAC2 o > Mr  IU  C7t?ll0  I L  AMM2  I  .  «*G?i  •  L  C 
11 

I  F  (  J CP  1  .  GT  •  ' f  P2  .  AM)  .  ‘.t p  1  •  GT  .  _  l  P.  |  u*  I  TE  C  7  * 2 1  2 2  )  .  f  ?  1 «  L  ■  C 1 
IF<  SEP1  -Li-,  tf  1  .Af..3.:tr  i.LT  .  GLr  )  Wr  I  T  L  (  /  •  21  20  :  l  PI  ,  L  Cl 
L  AM  til  J  -LAMBDA  <  1  » 

LAMd2l=LAMBDA«2> 

MAG 12=”AG11 
MAG22=*AG2l 
MAGI l  - wAGlll 
MAG  21  -  MAG  2!1 
St^2s5tPl 
SfPlaSFPO 

LCCl=CMPuA« f EQ.OElT A) 

2060  FC1MAT(/1X*«MAXI  HUM  VALUE  UtTkCTLO  FC?  L  AMdOA  C  1  1 1  •  /  3  X  «  •  t  AKUOA  Cl  >  = 
1  *«F12. 1 0  •  Jx»f  12. 1C  /JX ••MAGCLArnOAC i  )  i  -  •  ,  PI  £  •  « /<  X  .  •  *X  2  =  **Fd. 

2A/IKf«[)itrA  = 

2090  FOMHATC/l  X«  •MINIMUM  i/ALuE  DlTLCTO  FOf  L  AMHOA  C  1  )  :  * /3  X  ,  •  L  A"BD  AC  l  »  = 
1  •  •  F 1 2  •  1 0  •  3  X«F  1  2  *  1 )  /  3  X  «  •  *AG <  L *U  A  t  l  )  )  =  •  #H  0  .  •  /  1  X  »  •  \t  2  -  •«F*. 

24/5Xt»Gt.LfA  = 

2100  FORMAT  (/l  X**MAXlMliM  VALUE  DttLC't'D  F  OF  LAHBQAC  2  I  !  •  /3X*  *L  AMB0AC2  )  = 
1  ••Fl2«l0t3x«F12*10/.'X«*MAG(lA*'riOA(2ll  =  •  *F1C. -i/3  Xt*  VlG  = 

24/3  X,  •OUT  A  =  •*Fc../» 

2110  FUM  MATl/lX,  *  MINIMUM  V  ALLL  DETECTED  FCh  LAMBDA  C  2  )  !  •  /3X  ,  •  LAMBDA  C2  1  = 
1  •#F12#10#3X»F1A.10 /3X, -PAGCLA-hOAi:  )  )  =  •  *F  1 C  •  - /*  X  «  •  Kt  0  =  •«F.J. 

24/JX* -DELTA  =  *»F 

2120  FCPMAT  C/1X,*HCDl  St*°APA  1  JT.  PtAK  0£  TEC  T  fc  0:  « /JX  .  •  S  fPA  HAT  l  C!.  =  ••Flu 
l.**/ 3X  ••‘vfd  =  •»FwW'K»*UUlA  i  •  *  F  .:•*:/  I 

2130  FCr<  MA  f  C /l X  •  •  MODI  iFPAFA'IC*.  *IMNu*  Ut  rlCTFD:*/3Xt»  TfPA^.ATT'TN  =  • 
1F1J  .6/3X»»?.CQ  =  *  .Ft.  A/3X#«0t.uT  A  =  •*K6.b/) 

C 

C  PLOT  The  FIKST  gEVEn  EIGEWALU  • 

c 

20  70  NUMPL OT  =  7 

I  FC  M)r*G.LT,f.UMPLCTl  f  L  MF  LO  T  -  .‘.DEC* 

oo  Ji*;  xplctsi •mjmplct 

MAG  LAM  C  l  PLOT  I  =CAHf  C  L  AMii’JA  C I  r’LO '  )  ) 

I  FC  [t!K  >*E<3.  INFO  IX  Art  -  AFCPLC  T  >=*  k  Q 
IF  (  I6F  3*  k4E«  1NEQ )  XAn  °  A  YC  I  PL-  T  |sOtL  I A 

2140  COttUMut* 

C  Ai_ L  CU4WCCXAhPAY*MAGLAM*M  MFL  3  ?f-l  I 
OLL  T  A=  DEL  T  A+DELl  %*C 

2030  CONTINUE 

NEQ~MEG«uEQ!NC 

2020  CONTINUE 

GO  TO  2150 

77?  MRI  TEC7f2lF,u  ) 

2160  FO^MATC/l  XtlOO  (lh*  )  /  /lx••P^  uGnAM  WAS  T  Ef  MI*.  AT  ED  BECAUSE  NBIG  tfAS  C 
lUTiiUE  OF  T  TS  L  IMI  fi.»//t  X.  lOOCiH«M 

2130  CALL  ENOPLC17 
CALL  OwNEPL 
STuP 
EnO 
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APPENDIX  C 


Program  Listing  -  FOCAL 

The  code  FOCAL  performs  all  of  the  far  field  calcula¬ 
tions.  The  resonator  mode  is  propagated  to  the  far  field, 
the  centroid  is  located,  the  integrated  intensity  is  computed, 
and  the  beam  steering  angle  is  calculated.  The  output  includes 
intensity  profile  and  integrated  intensity  data  (TAPE 8 )  and 
plots  (TAPE 9) . 

The  input  data  is  read  from  TAPE7  and  is  entered 
interactively.  The  data  on  TAPE7  in  the  order  required  is: 

MAG:  cavity  magnification  (real) 

NEQ:  equivalent  Fresnel  number  (real) 

DELTA:  6  (real) 

MODE:  mode  number  (integer) 

ROOT:  mode  eigenvalue  (complex) 

XMIN :  minimum  x  value  at  which  the  resonator  mode 

mode  is  calculated  (real) 

XMAX:  maximum  x  value  at  which  the  resonator 

is  calculated  (real) 

INCX:  number  of  points  per  unit  x  at  which 

the  mode  is  calculated  (integer) 

NDATA:  total  number  of  points  at  which  the 

mode  is  calculated  (integer) 
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FIELD:  the  complex  field  values  (complex) 

SLOPE:  slope  of  the  phase  front  in  the  plane  of 

the  feedback  mirror  -  required  only  for 
mode  one  (real) 

Data  entered  interactively  includes: 

AMIN:  minimum  normalized  half  angle  value  at 

which  the  far  field  intensity  is  to  be 
calculated  (real) 

AMAX:  maximum  normalized  half  angle  value  at 

which  the  far  field  intensity  is  to  be 
calculated  (real) 

NUMF:  total  number  of  points  at  which  the  far 

field  intensity  is  to  be  calculated 
(integer) 

IRES:  response  to  a  posed  question  -  input  Y 

for  "yes"  and  any  other  letter  for  "no" 
CCPLOT56X  or  CCPLOT1038  is  required  to  generate  the 
Calcomp  plots. 


P  A  J  OR  AN  F  OC  AL  <  I  .PU!  .OUTPUT, 1  APt.  t  '.PUT  1 1  API  t  =  OU!PUJ  »  T  APT  T  ,  T  APt  ••  ,TA 
IP  E  J  I 


THIS  PROGRAM  TAKtS  T  Hf  NFA-t  FltLD  DATA  GENE'-ATEO  BY 
BAR  C 1  (FROM  IAPI/1  A  MJ  PHJPoGAIIS  Thl  h  OF  IS  THE 
FAR  FIELO.  A  SIHPi-L  IMEGfAl:..  ...uTlI.c,  OSES  THE 
FOURIER  TAAI.SFSIN  FES  THE  Pm  OPt.GA>  I  ON.  A  FLIT  EF  THE 
FAR  FIILO  llME'.SITY  PAITC«i.  t.  OINCRATtO.  1  Hf  Ilf  A N 
CENTROID  i*.  Inf.  LL-CATEU,  A'.U  I  HE  .  \  TEor  AU  G  I’.'TENSITt 
IS  CALCULATED  F.-OH  f HE  CENTrOID.  I  HE  If.TEGrATEO 
INTENSITY  IS  Al  .0  PLJTI£0. 

REQUIRED  II. PUTS: 


AH 1 1« 
AH  AX 

I 

cewra  l;pu  FCn  Per*: 
uppih  l  :hi  t  fcf  pl  . ’ s 

NUHF 

= 

TCTAL  •iUHIIi  '  OF  PC  I  f.  T  i»  AT 
IS  CAlCULAIi.0 

vhICh  the 

field 

IRt  3 

= 

RtcPCNU  TO  A  P ;> jL t)  aUE^TU*.  f*.Tff 
rts#  Ah  y  OTHtii  LETTER  F  Oh  KC). 

T  FOR 

THIS  C  JOE  Rl8Ui.;Ei  CCPLJIUJ-t  T  .A  CALC'. ME  ELCTS.  Oc 
CCPLOT-io*  MlK  Pr  £  V I  L  W I  NG  PLJTS  AT  A  GRAPHIC.' 

TERMINAL* 

FINAL  FORM:  J  (CT JhLS  IBS,?.  2LT  S.  M.  RJNALDI. 


COMMON  FIE  LOlHi  I  .X'lN  .XMAX.XI  .  C  .  CL  L  T  A,  ■  DA  T  A  ,  F  “  I  V  .  F  MA  X  ,  f  Tl C.NLMF.F 
IF  IN  r  loll  o>  ,f  EOC.M  VAo  uMi. .  u  )  ,r  I  THAI  l  Jc  >  .RAG 
COMPLEX  EYI  •F’.ILO,-'  iOT 
REAL  MAG  * NE  S 
DIMENSION  EAHELU7I 
OAT  A  LAHF.L/IHlOH  1/ 

EYE  =13. ,1 .1 

PI= J.lAlSOIcBJV  SR 

OATA  IYES/mY/ 

READ!  /.•IMAU.NEO.Ol  Ll  A.ROOl  <  «0U  T  ,  XM  !  N  ,  X«A  X,  'NCX.NOAl  A 
READ! ?»•> (FItLDlI )■  . = 1 t  ED  AT  A ) 

I  FI  MJOt.ta.l IF! AO  If .• I.LuPE 

CONDITION  FIELD  -  JET  PCNEK  :  1.  AND  BLANK  FIELD  OVER  FEEDBACK 
MIRROR. 

XINCM./I.NCX 
CALL  CONOI! 

OUTPUT  SUN  OATA  TO  CRT  ANO  cUTELT  FILE  ITAPf») 


as;  EE  <1.  .  1  C  I  HAG.I.CO,  til  L  r  A,x*t;  N, XMAX  .xIsC.NDATA  ,-C  >T 

UKI  ICIr.  ,10  >'-AG,  :  l  O,  DELI  A  ,  XN  I  >1,  A-<AX  .  X I NC  ,NOA  1 A  ,  Of  •  »  0J? 

10  F0RNAM/1X. •!-(..  PA-iAMtTlRS  (fil.s  DA, A  F  i  1 1.  »  *  // 1  X  •  *  <AG  -  .,F-.*/I 
1X..NE0  =  •  .1  •  .A/l  «.«0t  LTA  =  •  »F.  .A/XX.««M;>.  -  •  t  E  J  .«  /  IX  ,•  XMA  < 

2-  •  .FR.A/1X  ,.Xl>.C  =  •  .FR.A/lX,.  .MATA  =  •it]/lX,««  OE  »  •.I2/IX.* 

JMQUE  EIGENVALUE  :  •  »F  K  .2 • 2*  «F  1  j  .  /  I 

PRIM  MAXIMUM  ALLC  J  ABLE  SPATIAI  F  -  f  QUt  *  C  Y  5  “EA3  Of  ,  I  R  r  D  SPA  I  Al 
FREQUENCY  EIMITJ. 

All  H=  INCX/-  ..MAS 
MAI  ICI..JOM-AI  IM.AL IM 

•0  FOR  MAT!  IX  .«l  I  «l  l  S  ..  T  HI  Nu-MAl./LD  HAl  E  ANGLE  A-l:  «,F  .A,*  T.  •, 


I 
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2  joints* /ix  fit;  hi  Gt  a  rateo:*/ 1 

REAOC  S,*  )AMN,A*AX,  .UMF 
FMI  N=AMl*/C '  .*►  AG) 

FMAX-AHAX/ti  •  AC  I 
P  INC=  ||  h A  X-F  M  I  N ) /  (  •  ) 

CALL  FF CALC 

list  intensify*  if  utsi  to. 

WRI  T£  I  •»  *  A  3  I 

*0  F0-*  MAT  C /I  X*  *W„ClU  T  0  L  I  At  A  PtlMvUT  Af.O  C- T  LIST  I  *>  G  CF  THf  I’.'Ff. 

tsir y?»/i 
*EAO<S,lJo  >i-t 
\Q0  FvrtMAT(Al) 

IF!  IRE ;.Nt.  I Y r  SIGC  3  CC 
WRI  I£f»j**C  ) 

W*I  Tt  (- •  -i 0  > 

.0  FORMAT  </AX*  *\  HALF  A  \oL  l  •  *  I  X  t  •  i  M  F  I  TT  •  // > 

0  l  oO  1=1  •MM 

WRI  TE(o*70)FVALLL  (  l  >*FF  1ST  I  !  I 
WRI T£<J*70)FVAlLl C ),FF!NT(l > 

’3  FORMAr<‘jX.Fl»}.*jtbX*LlA. /) 

,0  COnTI SUE 

PLOT  INTENSITY 

j 0  CALL  PLQTStG.  *L  .,S) 

L ABEL ( 1  )  =  XChFA".  F1I.LU 
LAHfcH2)  =  K  hlf.d  r.i!  Iy 
LABEL! S)=10M  DILI  A  5 
LABEL*  1  J>=1  ;h  .»ni  = 

LAdtLUS)=10H  MuOt 

LABtL(t>l  =  UH  -AG  = 

LAHtLC  >I  =  1-;h  t.<  »■  f*.  HA 

LAHCLUOT  =1-JHLF  ANGLE 
LA3ELU1>=1CH  .CAi.rO  I 
LAUELU2)  MuMfiTl  SSI  TY 
ESCOOCClOt^tLAurLl  dlOELTA 
ENCODE  <10*1  Uu  *1  Artf  L  (HMMJ 
ESC  JPt  <10*110  «LA!UL  <lfc>  )MCUt 
ENCoOE<lO*liO*LABtL(b))MAG 
0  FORMAT  iF 4  •cj1 
100  FGhMATIM.*) 

110  FGRMAr<lXt»a*«lXvl?> 

120  FORMAT  <F4  •*  1 

call  HGKAPHCF VAlUL.FFIST *SUMF*L ABLL *1 *i*Q  > 

INITIALIZE  FOR  7Ht  ISTtGRATtO  IMk’.SITY  ROUTSE 
WRI  TE<*  *1 31  >Fl i C 

131  FGRMAM/l  K*«Pf  AK  FIELD  IS  LoCATtO  AT  A  =  •»F13.7/) 

WAX  TE  <  6  1 1 K  IF  t  SC 

130  FGRMAM/1X,*PI  AK  FIfLO  IS  LoCATtO  AT  A  =  •  *Flt’  .  7/ IX*  *Ff.  TC  R  MINIMUM 
l  AS  I)  MAXIMUM  HALF  A'  6lL  VALUIS  F  Oh  THt  l  N  T  l  G-  AT  f  0  •  / 1  X  *  •  I  S  T  F  %  S  I  T  Y  « 
20UTINC*  ANl)  ThI  SUM  HIM  -F  Pv.lM  ,  T  Z  BL  GF  M ATCo:  •/ ) 

RtAOC^**  >AM'S*A»AX*  .UMF 
WRI  TE  ( •*  *1  $2  )  AMi  AM  Ax  *f.LMF 

132  FORMATC/lX*-!*.!*  G-A '1  U  ISTl'.'ITY  Ri'JTlM  PA- A -f  T  f  *  S  :  •  /  3  X  «  •  A  "  I  *<  =  • 

ItFl  .5/  T  X«  «A*A  X  :  •  •  I  •«  .A/3X«  •  SLMF  -  ••13/> 

F  MI N  =  AMl  4/C. . «M AG > 

F  MA  X  =  A-AX/C.  .  *-AG) 

FISC=lf  MAX-tMl’.  I/I  IGHF-i.) 

CALL  FFCALC 


non  non  non 


CALL  P  [  mtltUOfCtMfcGI  > 

write  i.ntegfated  intensity  data  to  tape.>  and  output 

WR!  It  Ci,,l  NV  i 

150  FORMAT  C/1X,«U0  YOU  UA'iT  A  C*T  LlSTl?.*  .F  T  ML  INTEGRATED  I  f.Yl  NS  I  T  Y? 

!•/> 

READCS, J3JJ !PE$ 

Iff  ms •*«*..  :yes>go  to 
y  hi  I c c »i •  1  ;i  » 

170  FlH  MAT  (/I  X*  •  !Mi  Gf-  A  UD  »  N  It  *» :.!  i  Y  VALUE  ;.*/M  *•  f»  H  ANGLE  •  ,  1 1  X  *  •  I  NT  EG 
1RAT  tO*/  >K,«l  RLM  Ht  AM*»1  J*f*  I\Tt  r.ilTY  */c  X««ClNT»013*//» 

00  IrtO  1=1,1010 

u h i  ret'. ,i  >t  vai  ut  ( :>  ,p:  ntA< i  > 

190  F LIMA T <»X,E iO* 7*1 GX,F1j . H 
iho  cocriNut 
160  m*r  ret  « ,i  r«  ) 

00  200  1=1,1010 

yKI  rn  111  )  i  )F  VALUE  (  IlfPlTHAU) 

200  CONTINUE 

WRITE  INTEGRATED  1‘iTENSITY  U  AT  A  TO  THE  PLOT  FILE 

LAdtL<  i  >  =  iu  minti  gra  rto 

LArttL<2)  =  lCH  INTENSITY 
LABEL  Cl  1  >  =  UHlNTLGH  AT  LU 
LABEL (  1  2 1  =1 CH  INTENSITY 

CALL  MGkAPhIF VALUE .PITH A, IOIO*L  AHLl ,1  ,0  ,C > 

COMPUTE  l)IF  FRACTION  ANU  GEOMETRIC  lit  AM  EIEERING  ANGLES 

THETA  =  ULLTA*2*«MAG**.EQ*  CMAg-I*F/(NAC*1»> 

0MEGA  =  2.*Tf<F  ’’A/  (MAG-1  •  ) 

C£H  TROI  =  C  E  N  *  V  G I  -2  .♦  **AG 
«MT£4^,21(  ITht  IA,omf  GAiCF.  l  Kwl 
yR!  TEC.j  ,2 1C  >Th£I  A  »  J  Mt GA  ,CEV  MLl 

210  F0R"ATC/l  X,*KESL-NAT  cr  A\0  DIFFRACTION  ANGLE  SJ  •  /  3X  ,  •  GE  0"FT  A  I C  T!  LT 

1  ANGLE  =  •  ,F  1 A  •  7  ✓  3  X  *  ♦  N  L  i  0  N  A  T  J  R  Axis  TILI  =  •  »£  l  A  •  h/5Xt  •  01  FF*  ACT  I  ON 

2  ANGLE  =  •  tt  1A.  7) 

IF<  MOOt  .Ei3.1»«PIT€<^,22j  F^L.-PE 
ire  MuOE.EQ.  I  I  UR  I  TEC  ,22. ISL  .  PE 
220  FORMAT  COX  ,*  LOPE  =  *,E1A.7) 

CALL  PL  FTCCr.NNNM 
STOP 

end 

c 

. . . . . . . 

C 

SUHR  OUT INC  CONO! T 
C 

C  THIS  SUHRCUTIKE  CCNOITI.NS  T  Ml  fltLO  C  Wt  *  T  ME  PLANE  Of  THE 

C  FEEDBACK  MIRPUR.  SPECIFICALLY,  FIllO  VALUES  C\  T  Ht 

C  SURFACE  Akl  SET  TG  Cu„*..>.  TriL  INTEL  1 1 1 V  IS  THEN  NOFMALl/tO 

C  SO  THAT  IT-  INTEGRATED  VALUt  I-  l . : . 

C 

COMMON  F  ItLPCf.Cc  I  ,XMlf.,XMAX,Xl  .C  •  CELT  A  •  \OAT  A  ,F  M  I  V  •  f  m  AX  ,  T  I  NC,NU*F  ,F 
IF  IN  TIGS  Cl  #FL(C,F  VAlUE  (o  .  C  I  •  P  H  !t  A  C  i  S  ->  I  ,  *  AG 
COMPLEX  FILLO, tYE 
REAL  MAG 
P  1=3.1  MV*;  •  W.r 
EYE  =<0.,1 .1 
A2UPPEr =1  •♦DELTA 
A2L CM*  :-l.»Cf  LI  A 
X  =XHlN 
SUM -0. 

DO  10  I  =1  OA  T  A 
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i, 


nr»of>  on  ft  nnrinn 


I 


IF<  X.Uf  .  A  2i  .  yf*<  •  kU'J  •  X  .l  *  .A.  uKPl  r  >f  IElO<  il=«C.»rJ.) 

F*»AG2f  li.L0<  I  »  •CwNJG  (KIL'DUH 
IP= 11/21*2 

I M  [  *L  Q*  I  «i  -  •  l  « I  Q«!«iiATAIGC  *  t’  2\ 

IFC  fP.t'U.  1  »  UP=uU«*«».*F*"Afi 
XFC  IP.M.l>CUH-:UN»2«*F»AG 

GO  0  3U 

20  SlM-SUH*FMAG 

30  x=<»xr*c 

10  C*J  <  T  I.j'jF 

Sf: -UM*  *T NC/3. 

3Q-*  :Fs„Qi  T  l  f  > 

03  AO  I=1«M)ATA 
Fit  CDUUFULOd  j/:gk-f 
AO  CONTINUE 
RET  URN 
Z'.O 

c 

O*** ***• ••••• . . . . . . 

c 

subroutine  f  F  CALC 

rnU  OUUKCUTINL  PA  j  PAG  A 1  £  3  T  tt£  .UTPUT  F  .  LLO  TC  Irtt  FAs  FILLD, 

T  He  PC  AK  INTENSITY  POINT  I-  KETLhUEO  ALuNG  WITH  THl  FA-  F  I L  LO 
PATTErN. 

COMMON  F  I  EL  DC  ECU'  )  *X  mu,  t  XHAX  •  XI  ..C  .  Oc_  L  T  A  «  f.  DA  r  A  *F  N  I  N  *  f  H  AX  •  F  INC  .  VJHF  ,  F 
IF  I  '.  M b  J  0  >  *F  1 0 C  »  F  W  AL  UL  <  r,.’  e  >  •  »*  l  TH  A  C  1  3 1  )  .  P  AG 
C  OH  PL  t  X  FlUO.EYt  •jUH.l'.rAho,  I.PIF  ,:uhi 

A  CAl  hag 

F  -r  HP, 

HPl GHT -Z • 

E  Vt  M  j  .  .  1  .  > 

Pl=3. lAloVitSibc 

3ET  UP  00  LOOP  AOOUT  F  VALUES 

DO  10  IFUslfkUMF 
XMMIN 

I2PIF  =  — :VE*2.*PI-F 

SET  UP  OO  LOOP  T 0  T  AKE  FCUlILft  fHA\3FC.-lH. 

INTEGRA  T1  Cf.  TICmNIQUL  -  SlMPi,C\*3  NULE. 

SUH=C0.,0.> 

FVALUEC  (F0»=F*2.*MAG 
DO  20  l  X-  1  « f.OAT  A 
If.T  AftG-FICL  Dl IX )*CE  XPC lePir  *X> 

1  XP  =  UX/2)*2 

IF!  IX. FQ*  I.CH.lX.tO-MJATAIGu  T ?  3o 
IF(rXP.LQ.IX) SUH=SUM#a.*  IMAnG 
IFC IXP.NL  .1X1 SUH=SUH*2. •! M Ah  6 
GO  TO  A) 

30  SU1  =  $UHMM  Aft 6 
40  X  =  <  »Xl %C 
20  CCi  T I SUE 

S UH 1  =  SUH*  XI NC/3 • 

FFCMCIFO  l=  >UHl  «CCN JG(  :UH1) 

IFCFFIT.  rClFOI*LT»b^IbhT  I  GO  10  *i_ 

8P.IGHT-FFP.MIFC  ) 

FL0C=F*2. *MAG 
50  F=F  *FISC 
10  C  ON  T  t  N‘Jt 
RET  UK H 
END 
C 
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c**- . . . . . . . . . . . . . 

c 

SUBROUTINE  PITBIIDI  U*Cfr,TFu  l  I 
C 

C  T  Hi  S  SUBROUTINE  COMPUTE  :  iMtC-'ATF.J  INTENSITY  FC*  THE!  FAR  FIELD 

C  PAT  TEH'..  tit  A  H  C  E  M  l  L)  IS  T  AKt  N  Ai  THE  POINT  AT  UriICi  HALF  ThL 

C  TOTAL  POUEH  L  If  -  O'*  EITmEk  ^  £  01  - 

C 

COMMON  FltcUCtCC  I  tXMlNt>MAK  ,  XI  .C,  OlLT  A,  \0A  T A  . F p  I  S  •  F  «A X •  F  I  NC  •  NU“F  *F 
IF  IN  '  CijOCI  #FLCC»f  VAL  Uf  (*.w  J  )  *f  I  fMAC  I  5>  I  *  MAO 
COMPLEX  FIIU) 

UAL  HAG 
C 

c  utriHHiNt  total  pcgim  i^sior  limits  previously  (ntemeo. 
c 

SUM-0* 

or.  u  jm.slmf 
IP- i / 2 • * 

IFI  1  .EQ.l.wf  •  I«l  Q.NL'MF  )G0  T  2-.  * 

IFC  IP.e  U.  I  »  SUM- l UM* FT  I'.’Cll*4*. 

IFC  IP.;*t.I  loUM=SUM*f  FlN?  CU  *2. 

GC  10  1 0 

20  SUM  s  S5JM*FF  *  N  I C  I ) 

10  CClIlN'Jt: 

POytK=  ;ui«f  I  NC/ i. 

c 

C  LOCATE  tit-AM  CENTROID  WITH  AIO  wF  A  LXNlAR  INTERPOLATION# 

C 

POJtHA  =  PJWf  h*I.cv 
F  IyC2=F  I *«C •  I • 

F  IN  CA  -f  !  ’iC»  4  • 

SUM-C. 

I  -S 

MO  I=.  *1 

I  p  -  <  I  /:>>•? 

IFC  I  .CO. I  IGu  0  SO 
(Fcip.ro.  iiLUM=wUH*FFr.rc  1 1  *r  i.ca 
IFC  IP. NC.  1 1  SUM=f Uh»FF  In  T(  1 1  *f  tr,C2 
GO  TO  L3 

SO  SUH  =  jUH^FFINT  d  )«F  M  C 
60  IFC  3UM.LT.rCUf  RAIGO  >u  Nu 

If  <  I  P.t  (J.  I J  :,U*P  =  SUM-t  F;  \  T  Ci 

IFC  IP.M.I  1  SUMP  =  SUM -FF  IMC  I*F.  .Ci 

CESIHOI-FP.*  ♦ C 1 -2. )  *F I %C*  C  P  ‘A-lUMPI«F  INC/C  SU“-Iu  ‘P> 

Oli  T  -  CFMA  X-f  Hlf.)/Z  • 

FMAX  =  Ct./*Tr'r  I  ♦OlLT 
FM(  N  =  CFNT  *»!  t-DI'T 
FI.C-CFMAX-FM  INI/C‘iUHK-1  .  I 
CPiCf  n: h  H*  mag*. • 

URI  TC('..U,«  ICP 
WR;  TE  C  •  l  G T  2  )  CP 

1002  FCH*AM/1  X»*HEAM  CFMAu.D  L  C  A I  *  D  At  A  -  *.014.7) 

p o«*  t  h p  =  p  jy4  -  •  1C  : 

MR  I  TL  CN tl Ov 1 IPCGLFP 
yHI  Tf  C.>  .1  JU  l  IP'  #EHP 

1001  FO«MATC/lX,*Ptf LLNT AGE  *-GWl~  I*  REGION  USLU  FC*  LOCATING  CENTPOIO 
1=  •  i  F  .*»  •  S  I 
f  M*U«=?.*M4G*F:oC 
If RI  T£C  0  ,  \  U*F  J- 
U  H  l  TEC'  tloc  .H>rO- 

1000  F  G«  MAT  C /I  X» -MAX  EPRuH  I Cl\  ihu.O  LJCAIICN  CIS  VJ*M.  HALF  ANGLE  UN 
S(  Ti  >  =  ».t 14. 7) 

CALL  f  *  calc 

c 

C  INTEGRATE  ThL  INTENSITY 

r 


•  { 
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r»  o  r> 


i 


pit  ha 1 1 »=a. 

I  s* UMF  /2*  1 

FAC  nR-FINC/3. 

|N0t<=r/2 
I  di  u=im)«:x*i 
00  7  3  11-ltINOCX 
I A=  ?• 1 1 -2 
lB=r.  A*i 

ic=  ia*i 

PIT  BA(  II  ♦  1 1  =FACT  3K«  (FF  i  .T  U  ♦  l  Al*<i.»FFI  f.MI  ♦!«)  »FF  I  ,T  (  T  ♦IC  MTFIK’  (  I 
1- 1 A  1  •  FF  I M  (  I-  i  B  >  ♦  FF  I G  T  ( I  -  i  C 1  I  •  P  Z  T  H  A  (  1  1  1 

FVALULU1  )=<Xl-l.  >•»  UC*A.»MAC 

70  co.'iriGut 

fvalufunoi  x«f:*.c**.«mag 

ALT  IMN 

EGO 


10 


25 


27 

30 


C 

c 

c 


10 


20 


S  UK  K OUT  I  f*C  HOP AP MIX  ,  Y  •  N  •  I  D,  Go , f.P  ,  ;  1 

DIMCNSiGW  XdhYd)  ,i0(<?3)  t  IF  <%?.LQ.?>  CALL  PLCT  C-l  .  >5*2. 1C  .-3) 
IF  CNO.fQ.2)  GO  TO  Su  \  IF  (N).LT«.l  GO  TO  13 

CALL  SCALCi  X  ,7.  ,f««  1  1  V  CALL  ,C  ALL  (  Y  ,3.  ,N  ,1  1 

CALL  PLOT  (l.,ll.,2l  t  CALL  Pi_  L  T  l  <•  .  5  ,  U  .  ,  2  ) 

CALL  PLOY (M  .5  *u  .  ,  ?  1  %  CALL  PL 0 7 ( J . , 0 . ,2 1 

CALL  PLOT  (1.35,  1.3*>,-31  1  CALL  PL  GT  C  - .  .30  ,-2  > 

IF<  IOU  1.CQ.9V?  1  GO  TO  c5 

CALC  PLOT (.1 ,-.1,-3 1  S  CALL  PL C T ( 0 •  . -2 • « *2 1 

CALL  SYMBOL < .25#. 3* .0  7 ,  1  D(  1  1  ,  V  •  ,T  J  1 
CALL  SYMBOL  (.45, .3,  .0  /»  0 « 3  >  ,Vi  .  ,1  J  ) 

CALL  $YM»J0L(.t5,.3,.v  7,0(131,  -  •  ,2  j  ) 

CALL  SYMBOL  C.f 5,.3* .C /. ! D< 1  I*  . *2  > 

CALL  SYMtTCL<i.05*.3  f.C /*ID(‘- I,  > 

CALL  SYMBOL  (1.1  5*.  3  ,.G  7  ,1U(  /  >,  >„  .  .201 

CALL  PLUT<C  «'C.  ,31  »  CAlL  PLc  T  <  1 . 25  *0  •  ,  2  > 

CALL  PLOT (1.25,2. ,2  1  \  CALL  PLC T ( 0 . ,2 . , - 2 1 

CALL  PLOT (-.!,. 1,-51 
CALL  PLOT <5.4 1 

CALL  PLoT<u.,-:.3G,-2)  J  CALL  PLOT  (-3.  stC  •  ,-2  1 
CALL  PLOr ( 5.3 ,. 75,-  M 

CALL  AXI5(0.,0.,ID('*),-J0,7.,'/0.,X(N^1),X(T. ♦211 
I F( \3.lQ. 31  GC  TO  2 / 

CALL  AXIS(0.,G.  ,  I  D  (  11  >,.C,5.  ,1?  •  »  Y  (  ?«♦  1  1,Y(L*21  1 

GO  TO  33 

call  axis(g.,o.,io(wi,:'g,3.,i.  .  ,y cn*  i > ,t <»,♦?>  > 

Y  Cj*2>--Y  (WJ  t  CALL  L  Z  GE  (  Y  ,  X  , '« •  1  ,  *»  P  ,  *1 "  1 

Y(S *2  l=-Y (0*21  t  CALL  PLL T ( I.c5,-2.1C,-3| 

ACTUKS  »  LgO 


SUTcOUTIgE  AXICIXC, Y:  ,L,GC»  «L.A  .G  ,  P  “  I  0  •  DM 

01  MEGS  1  OS  LU1  1  A  -  A  \  G  •  .5  •  l  *  1 5  *»  S  .  u  •  *  Dx s . 1 • CO. ( A  1  l  OY  =  .  1  I M  A 1 

r  c=i  Si’G  ,  vc )  i  ,c  =  iA»ii  ,c  >  i  >-=.i  »  .=i  *  x  =  x-  i  y-y:» 

CALL  PL  jT  (  X  ,  Y  ,3  1  1  X-X^UX  l  T  =  f*Jf  *  CALL  PLCT(X,Y,/> 

CALL  PLOT (X-.21 •  0  Y*  iC,Y*.21  *JX«  C,J  » 

I  F<  l.f  •).?.)  CALL  PL1’  (X-.«2«UY*iC,  Y*.A2*|IX*:C,21 
I  F(  *i.L  0 . 1  0  I  CALL  PL  *'  Y  C  X  -  •  /  ;  •  OY  •  »  C  , Y  I5»DX*IC,<.) 

N -M  ")0  (  '<  ,  l  j  1  ♦  1  S  1  *.1  <  H.L  '  •  -<l  >  G  TO  10 

AsANG*(  i  C  ♦  1  1*45.  1  CX=lo.*iU  %  DY  =  U.*UY 
C  =-  .175».125*XC  1  U  =  .l-*.3i,«IC 

x=x : ♦c*ux-u*dy  i  y • r . *c*ur jx 

R  =A  MAX  1  (  AOj  (H HIM  ,  A iiL  (  «  L  1  1  1  k - AL  GO  1 0 ( J  1 

XRs  IGMAHSO  >  1  1  IFCr.L  ...  1  i  ---«!*-♦  1  I  1  I  -  =  I  <  -  M  ^0  (  I  a  ,  ?  1 

RI:*M1*«/U.*«  IH  S  11  1  =U- /1C  .••;»  I  =L. 

EGCOOE* /,  1C  1  ,  Slrl  S  CALL  SY  M.ICl  4  X  ,Y  ,  .C  7  ,  i  ,  A  •  7  )  1  ^lsPl*0Jl 
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X  =  X*i)«  t  T-Y*nr  %  1  ♦(-.Lt.-LI  «.J  20 

it  =C-L-.  1*V  C|/4  .  %  (.5.1*.  -1C 

x=xo*K.ox~c*ur  t  Y-r  •  ►•oy*c«i:x 

CALL  SfMII  H  (X,T,.l  ,L,A  V.,*.  tl  *  IFO.tU.  I  U  ▼ 

Cf.Ci’OC  C  -»•  |.  .**  :  »  STALL  T"t«  .LC  •••.It*  ,Af.G«N) 

call  y^U(X|YtA) 

trie  JOt  (  3  •  1 L  5«  SI  I*  S  CALL  -  THoU  (X  »  Y»  .  :  /  9S*A\r>  ,  3  ) 

101 

102  F^-AJhH  *13  » 

103  FC*mAT<!3» 

HCTUkN  \  i  .!) 

C 

C  GUHrtft,UL[  <C  ALII  DATA  .LENGTH, T.,K  > 

c 

C  Kf  AL  DATA  =  K.t-  L' !  Ft  A.  MO  A-hAY  Lf  UA  T  A  T  *  H:  ;CALfO 

C  ;?jTtbt*  t  i  t.  ufm  f  '  »  DA  r  A  b  . ;  . ;  ; 

C  ‘-LAl  Ltt  GtH  =  L*‘gTh  J  •  HL  rLv!  Axl.  (t  iOt  .  *.  l-.CHL ; ) 

C  INTEGER  K  -  UMAU  lA«A‘i’t«  i.CLUDlu  ft-  Cl  *•  P  A  *  !  H  !  L  .  T  Y 

C  WITH  rut  £  QO  1  V  AL  r  *•  T  LALC'MP  ;.lihR  CUT  I*.C 

c 

C  Trie  F  o  LL  >-  W  I  NO  VALUL  A-E  i>£!oM,IU: 

c 

C  OA1ACF.M  I  3  AUJu.Uu  DATA  MIM^U* 

C  IIAVA(N*2)  =  "MU*  .CALL  ULTw-.  !\  UA  T  A  UMTS 

C  Ft*  lENGlM  U.IT  (E.G.  V^LT^/I’.Crtl 

c 

. . . . . . . . . 

c 

SUir  JUflNfc  '  CAU  t  OA  ’  A  ,Li  .0  I  h  v  K  > 

■UAL  OATA('),  it  NOT  r*  •  .1  <‘,| 

OAT  A  A  /l.  •  2.  *  .  i  \. .  1  .  •  / 

C  C(.MPJTl  TMt  nAU  SCALE  FACt  - 

□  Hi  '.SOHAX  =  OA  f  A  (  1  I 

o:  u  i=i.. 

!F( DAT  A ( !  >  *L  '«  J-IM  UM.  .  =  OATAU) 

IF  (OATA(I).  Gl  .  u“A* )  O’*  A  x  =  DA!  4(  :  I 
13  Continue 

C 

C  EXCLUDE  HI  VIAL  Lh4(.*  C A  .1  G 

C 

OAT  A  <  N  ♦  1  I  =  OF*  I*. 

0  A’  A  C  '*♦  2  I  =  1.0 

IF  (LC  %GT  H  .LC.  J  .J  .0...  UMAX  .L  Q.  OKlf.  >  tfTU»-N 
RAJ0F  =  ( UP  AX  -  OMIM  /  LENGTH 

C  R  Ay  jF  =  3  F  KANT  •  10.  ••  SFCXP*  nMt  -  L  1  .LC.  OFPAf.T  .LT.  1C 

SFt  XP  =  AIM!  AL06 1  i(  HAWSE  )  I 

IF  (  R  A  V  G  F  .LT.  l.J  )  ^FLXF  =  i  L  XP  -  1.3 

SFMANT  =  RA«3F  •  10  .1  ••  (*vFtxH 

C  LOCATE  NEXT  LARGER  "MCE*  3CALE  FAC'CR 

O  G  2  U  I  =  1  «  -> 

20  IF  (  GF(I»  .GT.  SFMAI.T  I  Gv,  Tc  3c 

PPI\T*«*  *»C  AL  £  :  SCALf  FACT,*  Er->  ...  ■  i  *fTU*N 
30  SfMCE  =  .rCl)  •  1J  .  ••  oF:<P 

C  CJMPJTE  A  □  JUS  TE  O  DATA  MIMMLH 

aujmi\  ^  a:.,t  (  dp;*.  /  fmcl  t  • 

if  (  Af)jHi,.  .or.  um:n  i  adj«:i.  -  aojmil  -  •4fmc‘ 

IF  (  (UMAX  -  AUJMI*.  I  /  -FMCE  .LE.  LENGTH  >  G 1  T.  A 2 


c 

Nff  J 

r ;  u:t  thf 

\LXT  UA«GL-  jCAlL 

F  ACT  o.^ 

IF  < 

;  .L7.  5 

1 

sf  m  ci*  =  ;fii*d 

•  1  0  •  j  •  • 

CFLXP 

IF  ( 

I  .to.  5 

> 

SFMCt  =  2..: 

•  10  mi'  •• 

ZFC*f> 

ADJOIN  =  aim 

( 

OF I  \  /  IMCt  )  • 

Met 

if  c  aojm:*.  .gt.  on:m  aojn;,.  =  adj.^n  -  \fmce 

40  cc.Tiv  r. 

0Ar  A(  «»l I  -  ADJFlf. 

OATAIN*.*)  =  JMCt 

Af  f  UAfc 

E'.J 
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